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LIST OF SYMBOLS, ACRONYMS AND ABBREVIATION 


name of the coordinate system, B denotes body-axis, W denotes wind-axis, A 
denotes air-mass axis, and U denotes inertial-axis 
rotation matrix from coordinate system {A} to coordinate svstem {C} 


skew symmetric matrix 

forces in {C} 

moments in {C} 

position of point D, measured in {C} and expressed in {C} 
position of point D, measured in {C} and expressed in {A} where 
“GE =¢ IIe )z°P, 


‘position error 


velocity of point D, measured in {C} and expressed in {C} 
angular velocity of point D, measured in {C} and expressed in {C} 


time derivatives in the body-axis {B} 


time derivatives in the inertial-axis {U} 
distance 

length 

atmospheric density 

dynamic pressure 


weight of i 

mass 

mass matnx 

apparent mass coefficient 

apparent mass matnx 

moment of inertia 

drag 

reference area 

drag area 

roll 

pitch 

vaw 

roll rate 

pitch rate 

yaw rate 

glide ratio 

control deflection 

forces in the x-, v-, Z-axis, respectively 
moments in the x-, v-, Z-axis, respectively 
change in force with respect to 7 

change in moment with respect to 7 
Capital letters/symbols denote vectors or matrices 
Small letters/svmbols denote scalars 

If coordinate axis symbol is omitted, {U} is assumed 
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Subscripts: 


cae center of gravity 

| suspension lines 

C canopy 

p payload 

a actuator 

SS steady state 

Xx or | X-ax1S 

yy or 2 y-ax1S 

zz or 3 Z-ax1S 

e error 

AERO aerodynamic 

AM apparent mass 

GRAV gravity 

Abbreviations: 

A-D Analog to Digital 

AGAS Affordable Guided Air Delivery System 
AHRS Attitude Heading Reference System 
C/A-Code Course Acquisition Code 

CARP Computed Air Release Point 

GER Circular Error Probable 

DOF Degree of Freedom 

GNC Guidance, Navigation, and Control 
P-code Precision Code 

PMA Pneumatic Muscle Actuator 

PPS Pulse Per Second 

RAWIN Radiosonde Wind Measuring System 
RMS Root Mean Square 
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| INTRODUCTION 


The United States Air Force Science Advisory board was tasked to develop a 
forecast of the requirements for the most advanced air and space ideas to project the Air 
Force into the next century. The study, encompassing all aspects of Air Force operations, 
assessed a variety of technology developments critical to the Air Force mission. This 
study culminated in a report titled "New World Vistas, Air and Space Power for the 21st 
Century."’ The study identified a critical need to improve the Point-of-Use Delivery; that 
is, getting the materiel where it needs to be, when it needs to be there. Airdrop is an 
important aspect of Point-of-Use Delivery. The report indicated that immediate 
ete vements are needed with emphasis provided by the statement: "In the future, the 
problem of airdrop should be treated as seriously as the problem of bomb drop." 

To date, significant emphasis has been placed on the development of large-scale 
parafoil systems. These systems provide the accuracy required with delivery from high 
altitude and large offset distances. The drawback is prohibitive cost for each pound of 
payload delivered. Alternate approaches were required to reduce system cost. The team 
of the United States Army and Air Force, The Boeing Company, and Vertigo, 
Incorporated is evaluating alternative airdrop technologies. These efforts include the 
design and development of the Affordable Guided Airdrop System, which incorporates a 
low-cost guidance, navigation, and control system into fielded cargo air delivery systems. 
This study focused on evaluating the feasibility. of the AGAS concept and encompassed 
the design and execution of a flight test program to assess dynamic response of a flat 
circular parachute, to the design of initial guidance and control techniques, and to 


evaluate the feasibility of the AGAS concept. 





i. SYSTEM DESCRIPTION 


A. AFFORDABLE GUIDED AIRDROP SYSTEM (AGAS) 


The Affordable Guided Airdrop System? (AGAS) is being evaluated as a low-cost 
alternative for meeting the military's requirements for precision airdrop. Designed to 
bridge the gap between expensive high glide ratio parafoil systems and uncontrolled 
(ballistic) round parachutes, the AGAS concept offers the benefits of high altitude 
parachute releases but cannot provide the same level of offset from the desired impact 
point (IP) as high-glide systems. The design goal of the AGAS development is to 
provide a Guidance, Navigation, and Control (GNC) system that can be placed in-line 
with existing fielded cargo parachute systems (G-12 and G-11) and standard delivery 
containers (A-22). The system is required to provide an accuracy of 100 meters, Circular 
Error Probable (CEP), with a desired goal of SO meters CEP. No changes to the 
parachute or cargo system are allowed. 

The current design concept includes implementation of a commercial Global 
Positioning System (GPS) receiver and a heading reference as the navigation sensors, a 
guidance computer to determine and activate the desired control input, and the 
application of Pneumatic Muscle Actuators (PMAs) to effect the control. The navigation 
system and guidance computer would be secured to the existing container delivery 
system while the PMAs would be attached to each of four parachute risers and to the 
container. Figure 1 illustrates the concept. Control 1s affected by lengthening a single or 
two adjacent actuators. The parachute deforms creating an asymmetrical shape, 


essentially shifting the center of pressure, and providing a drive or slip condition. Upon 
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deployment of the system from the aircraft, the guidance computer would steer the 
system to a pre-planned trajectory. This concept relies on the ability of sufficient drive to 
be produced to overcome errors in wind estimation and the point of release of the system 


from the aircraft. 


GNC package 





Figure 1. Affordable Guided Airdrop System” 
B. PARACHUTE 
The C-9 parachute was selected for this feasibility demonstration due to its 
availability and representation of the larger cargo-type parachutes (G-11 and G-12) on 
which this system will ultimately be used. Although the C-9 was initially designed as an 


ejection seat parachute, it is a standard flat circular parachute as are the larger G-11 and 


G-12 cargo parachutes. A flat circular parachute is one that when laid out on the ground 


forms a circle. 





Figure 2. C-9 Parachute 


. 


The reference diameters of these chutes are 28 feet (C-9), 64 feet (G-12), and 100 
feet (G-11). The reference area of the C-9 parachute is taken to be the total surface area 
of the canopy (a circle of 28 foot diameter) and is 615.8 square feet. The C-9 is static- 
line deployed and utilizes 28 suspension lines connecting to four risers. 

A cargo box was suspended from the system and housed the remote control 


system, control actuators, and instrumentation system. 


Ge ACTUATORS 


Vertigo, Incorporated developed Pneumatic Muscle Actuators’ (PMAs) to effect 
the control inputs for this system. The PMAs are braided fiber tubes with neoprene inner 
sleeves that can be pressurized. Uninflated PMAs as installed on a scaled system are 
shown in Figure 3. Upon pressurization, the PMAs contract in length and expand in 


diameter. 





With four independently controlled actuators, two of which can be activated 
simultaneously, eight different control inputs can be affected. For this demonstration, a 
throw of approximately 3 feet was selected [measured to be 3.25 feet during ground 


testing]. When depressurized, the PMAs are completely flexible allowing for efficient 


packing of the actuators with the parachute. 





Figure 4. Packing the Parachute and Actuators 


The concept employed for the AGAS is to fully pressurize all actuators upon 
successful deployment of the parachute. To affect control of the system, one or two 
actuators are depressurized. This action "deforms" the parachute creating drive in the 
opposite direction of the control action. Figure 5 illustrates this action while the Figure 6 


illustrates the parachute deformation upon control actuation. 
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Figure 6. Parachute Response to Control Input (video clip) 


D. MASS AND CENTER OF GRAVITY 


The weights of the major system components were determined through direct 
measurement or calculation using known material weights (e.g. parachute canopy fabric 
weight of 1.1 ounces per square yard). The major system component weights are 
summarized in Table 1. 

The theoretical center of gravity of the system was found to be 11.4 feet vertically 
upward from the center of the payload for a standard sea-level atmosphere. This 
calculation includes the mass of the trapped air in the canopy and as such is dependent on 


the air density. The relationship of center of gravity to air density is as follows: 


(1) 1, =[1,(oVg+W,)+1W, +1W, +1LW,\VW 


System 


Substituting known parameters equation 1 reduces to: 


0 0 
1, = 4323p +0.5 feet and P, =| 0 |= 0 
l| [43230 +.5 


Weight (pounds) 
Component sea level - standard atmosphere 
Parachute 


Actuators (total for four) 


Payload 
cargo box 
instrumentation 
actuator valves/nitrogen tanks 


Total System 


Table 1. System Weight 





E. MOMENTS OF INERTIA 


The moments of inertia of the parachute (including trapped air mass), payload, 


suspension lines, and actuators were determined as shown below. A summary of these 


values 1s shown in Table 2. Due to symmetry of the parachute and payload system, all 
cross products of inertia where assumed to be zero. 


Moment of Inertia (slugs: ft’) 


Component so —_ - standard atmosphere 
a 177 | 
|SuspensionLines neglected 
negiected 


Payload 5.) 328. 39.9 
cargo box 


instrumentation 
actuator valves/nitrogen tanks 


Total System | _3194.2| 3194.2 


Table 2. Moments of Inertia 





i Parachute 


When inflated, a flat circular parachute approximates a hemisphere with a radius 


(r). The moments of inertia can be found as follows: 


5 
“= 


(2) [ye ely Hl = mr , 


These values then must be translated to the principal axis located at the center of 


gravity at the system. The parallel axis theorem” is applied as follows: 


Cr 1, nd, and eerie 


i 


(4) I RS > +4 md,” where: d, =d, —1,,. 


XX YY 


eS) 


The moment of inertia, like the center of gravity, varies with air density as the 
trapped air mass‘1s a significant contributer to this term. Substituting m=m, + plVg 
(where m, is the mass of the parachute material and all known parameters: 

Tye = Typ = 7344369 + 1449.0 


I, = 831529 + 53.5 


2 Payload 

The cargo box dimensions were 2 feet by 2 feet by 2 feet. It is assumed that the 
center of gravity of this box was located at the center of the box. The moments of inertia 
can be found as follows: 
(Smee = 1 ae = <m(a* +/°); where a=. 


Again applying the parallel axis theorem, the moments of inertia about the 


principal axis were found. The parallel axis theorem 1s applied as follows: 


(1 =) a1 pe dagen ie 


3; Actuators and Suspension Lines 


The actuators and suspension lines were treated as slender rods with the body axis 


located at the end of the rod. The moments of inertia were found as follows: 
(7) Lge = Ty = = and Ip — 0. 


These values for this system were found to be extremely small as compared to the 


values for the parachute and payload and therefore were considered negligible. 
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Ii. CONTROL SYSTEM 


For an airdrop mission, the aircrew will determine the Computed Air Release 
Point (CARP) based on the best wind estimate available at that time. The aircraft will 
then be navigated to that point for air delivery of the materiel. Should the wind estimate 
and calculation of the predicted release point be perfect and the aircrew gets the aircraft 
to the precise release point, then the parachute would fly precisely to the target without 
control inputs: However, wind estimation is far from a precise science. The calculation 
of the CARP relies on less than perfect estimates of parachute aerodynamics and the 
flight crews cannot possibly precisely hit the predicted release point for each airdrop 
mission. Therefore, the AGAS control system design must help overcome these potential 
errors. 

A glide ratio (ratio of horizontal to vertical velocity) demonstrated in flight test 
was approximately 0.4-0.5. Considering this relatively low glide ratio and a descent rate 
of approximately 25 feet per second, the AGAS can overcome only a twelve foot per 
second (approximately 7 knots) wind. It 1s therefore imperative to implement the system 
to overcome poor estimates in the wind and not try to steer the system against the entire 
wind. In other words, the drive of the system is insufficient to attempt to fly straight to 
the target but is sufficient to overcome up to a twelve foot per second error in the wind 
estimate. For this reason, a trajectory control approach was selected. 

A pre-planned trajectory, based on the best wind estimate available, must be 
determined and provided to the guidance computer. The GPS navigation system will 


provide a continuous position of the system. The guidance computer will compare the 


actual horizontal position, at the system's current altitude, to the planned trajectory at that 
altitude. This represents the position error (P,) at the current time. 

A tolerance cone is established about the planned trajectory (Figure 7) starting at 
600 feet at the beginning of the trajectory and gradually decreasing to 60 feet at ground 
level. Should the position error be outside this tolerance, a control 1s activated to slip the 
system back to the planned trajectory. When the system-is within 30 feet of the planned 
trajectory the control is disabled and the parachute drifts with the wind. Thirty feet was 


selected to encompass approximately 1-sigma of the GPS errors (Selective Availability 


off). 


no control 


¢ Generate predicted trajectory based on wind estimate 
¢ Establish control tolerance 

¢ Compare actual position to predicted trajectory 

* ‘Drive’ to predicated trajectory 





Figure 7. Control Concept (double click to see animation) 


As outlined above, the control system relies on the current horizontal position 
error to determine if control input is required. This position error (P.) is determined in 
inertial space and is then rotated to the body axis using an Euler angle rotation with 
heading only. The resultant body-axis error (Pp) is then used to identify which control 


input must be activated as shown in Equation 8. Two components are returned, a + or - 


2 


for the x-axis and a+ or - for the y-axis. It was assumed for this simulation that a +x 
would activate control A, a -x activates control C, a +y activates control B, while a -y 
activates control D. The actual rigging of the operational system must align these control 
actuators to the compass reference line to ensure proper control. Control A is assumed to 


be aligned with the compass zero reference line. 


; le 
(8) input = sign IP] 
b 


The magnitude of this calculation is used to determine if the selected control input 
will be activated. If the maciitide is greater than 0.3, then that control is activated. This 
will allow the activation of a single control input or two simultaneous control inputs. 
Preventing activation of control when this value is less than 0.3 will ensure that 
unnecessary control inputs are not activated when the predominate error is in a single 


direction. Figure 8 illustrates the region of active control. 


compass zero 


reference <a 
controlD | > 


activated direction of 
predominate 
error 





Figure 8. Control Activation 
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IV. TESTING 


The objectives of the test program were to investigate the apparent mass effects 
on flat circular parachutes and to obtain sufficient data to assess the feasibility of the 
AGAS concept. To meet these objectives, ground testing of the C-9 parachute and flight 
testing of the prototype AGAS, using remotely controlled activation of the actuators, 


were conducted. 


A. INVESTIGATION OF APPARENT MASS EFFECTS 


To illustrate the effects of apparent mass on the dynamics of a parachute, a test 
program was established to collect the forces along the z-axis of the parachute during 
acceleration. The concept defined by Vertigo, Incorporated” to evaluate the apparent 
mass terms was implemented. The C-9 parachute was attached to a tower installed on the 


vehicle with two attaching risers and towed behind a ground vehicle (Figure 9). 





\ Gs, 


The vehicle first accelerated to a speed of approximately 10 feet per second. 
After a brief period of constant velocity, the vehicle then accelerated at a near constant 


acceleration. Figure 10 shows a typical velocity profile for these tests. 
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Figure 10. Tow Test Velocity Profile 


The forces in the attaching risers were measured with strain gauges. In addition, a 
"box" was constructed on top of the tower to measure the vertical and lateral forces on 
the riser attachment points. By measuring these forces, an estimate of the parachute 
angle of incidence could be derived. Figure 11 illustrates the instrumentation 
configuration at one riser attachment point. The velocity was measured with differential 
GPS and an anemometer. 

To estimate the actual apparent mass parameters, the motion of the parachute 


would be ideally fixed in one direction. In the Vertigo test, the parachute used was very 
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stable. Therefore, the assumption of motion in only one axis was very reasonable. With 


the C-9 parachute significant oscillations were observed. 





Figure 11. Instrumentation Configuration: Riser Attachment 
A plot of the calculated incidence angle (0) demonstrates the oscillatory nature of 
this parachute. Therefore, estimation of the apparent mass coefficients using these 


techniques was not possible. 
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Figure 12. Incidence Angle (0) 
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Even though precise apparent mass coefficients cannot be found, insight to the 
apparent mass contributions can be gained by analyzing these data. Figure 13 
demonstrates the influence of apparent mass on parachute dynamics. The total force 
along the parachutes z-axis, as measured in the risers, 1s plotted along with the steady 


state drag calculations as found by: 


QO D,= = pV -S-C,; where So is the reference area of the parachute and Cp is 


the steady-state drag coefficient (0.68). 
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Figure 13. Apparent Mass Effects on Parachute Dynamics 


The figure shows increased measured force above the predicted steady-state drag 
along the parachute when the parachute is Accelerate! Note during the first 
accelerations that the measured force did not rise above the steady state drag estimate. 
This is likely due to the parachute inflation and the greatly changing angle of incidence 


during this inflation. Concentrating on the area where the parachute has settled into pure 


18 


oscillations, the apparent mass effects can be observed. Figure 14 presents an excerpt of 


data for this time: 


——— steady state force 
seer measured force 


wa 
% 
sr) 
en oe 
=} 
an 
as 
a 
Cs 
a 
Or 
i> 


ee, Time isecdads). fe 





Figure 14. Apparent Mass Effects 


These data demonstrate that the apparent mass effects must be considered in flight 
dynamic modeling. Even with the relatively small accelerations experienced here, the 
effect is very large and varies greatly with acceleration of the parachute (as observed with 
the oscillating parachute). Techniques for estimating apparent mass coefficients are 


presented in the system modeling section of this report. 


B. FLIGHT TESTING 


The flight test effort focused on the collection of flight dynamic data to support 
modeling of the AGAS concept. The flight test effort was conducted with four actuators 
in-line with a C-9 parachute and a one-half scale container delivery system. The 


actuators were activated using a manual radio control system. Flight dynamic data were 
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obtained including the position, velocity, acceleration, attitude, and attitude rates of the 
system. It was necessary to correlate these data with control inputs. Therefore, the state 
of control activation was monitored. Parachute performance is significantly influenced 


by the winds. It was critical to this effort to measure the winds as precisely as possible. 


C. INSTRUMENTATION 


Ideally, both the parachute and payload would have been instrumented to collect 
all necessary data. However, the state-of-the-art in instrumentation 1s not yet sufficient to 
adequately instrument the parachute itself. As a result, only the payload could be 
instrumented. A custom instrumentation system was developed and it included a 
differential GPS system for precise position and velocity, 3-axis accelerometers for 
acceleration, and an Attitude Heading Reference System (AHRS) for 3-axis attitudes and 
attitude rates. Pressure transducers were put in line with the pneumatic actuators to 
monitor their action. Figure 15 illustrates the instrumentation design. A summary of the 


major instrumentation components is provided in Appendix D. 










Sensor Instrumentation GPS Instrumentation 


Accelerometer Recorder 
Tnad 
Analog to Processor Processor 
Digital CPU1 CPU2 
Pressure Recorder GPS Receiver 
Transducers 


Figure 15. Instrumentation Block Diagram 
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The critical part of this design is the time synchronization of data from all 
sources. To achieve this, the IRIG time generator was synchronized with GPS time using 
the Havequick Time Interface and the 1-Pulse Per Second time sync. The AHRS data 
were available at its interface approximately 71 times per second. This is the fastest rate 
for all the sensors. Therefore, these data were used as the key for capturing data from all 
sensors. When an AHRS message was first detected by CPU1, the first record was time- 
tagged and the entire message was sent to the recorder. At that time, the Analog-Digital 
(A-D) converter was polled for its data. The first record from the A-D converter was 
time-tagged and the message sent to the data recorder. The GPS data from the 
differential receiver were independently captured (with an embedded time-tagging using 
GPS time) and recorded on a separate PC-Card. The synchronization of the data was 
validated in two ways. First, it was desired to have a discrete event that would effect all 
data sensors. The obvious event was ground impact, which resulted in immediate 
changes in the data from each sensor. This discrete event showed that the data from the 
independent sensors were time synchronized to less than 200 milliseconds (the rate of 
GPS data). Next, the acceleration data were integrated to estimate velocity. This 
estimated velocity was then compared to the GPS velocity (after rotation from the body 
reference frame to the inertial reference frame using the AHRS attitude data). These 


results also showed the data were time-tagged to within 200 milliseconds. 


D. WIND ESTIMATION 
Two methods of wind estimation are presented. The first includes the accepted 
standard of the Radiosonde Wind Measuring System (RAWIN) system used throughout 


the test community. A RAWIN balloon was launched at approximately one-hour 
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intervals near the release time in the vicinity of the Drop Zone. This system provided a 
direct measurement of the winds as a function of altitude. Although an accepted 
standard, the RAWIN system has limitations in airdrop operations. The largest problem 
is that real-time winds are not available. The balloon must be launched and data 
processed resulting in approximately a one-hour delay. Figure 16 illustrates the 
magnitude of wind changes over time. Data from three RAWIN launches are presented 
as a function of altitude. It is clear that as the parachute gets close to the ground, the 


wind changes can be significant. 
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Figure 16. Wind Changes Over Time 


A second method of deploying a "calibration system" just prior to release of the 


test payload was implemented. A tri-lobe parachute was used with a reference drag area 


Jed 


(CpS) of 90.88. Oscillations of this parachute were observed to be very small. The 
system was weighted at 58 pounds to provide approximately the same descent rate as the 
C-9 system. The calibration system was instrumented with differential GPS on the 
payload. Initially, the wind estimate was simply taken as the ground velocity as 
measured by the GPS. This approach effectively considers the calibration system 
massless and therefore does not account for changes in momentum. 

To validate the ability to simply use the measured ground track velocity as the 
wind ee a model of the calibration parachute system was developed. A point-mass 
system was assumed with the only forces on the system being drag and weight. The 


applicable equations of motion are: 


(10) X =(m+a@,,)u =—Deosy cosy 
(11) Y=(m+a,,)v =—Dcosy siny 
(12) Z=(m+a,,)w=—-Dsiny +W 
where: W is the calibration system weight, D is drag (D=qC,S), vy and y are 


the flight path angle and yaw angle respectfully. The reference angles are: 


(13)  siny=——-; cosy = a... siny = “cosy = —— 
a WEE oe 


Substituting, rearranging terms, and putting in state space form: 


S 


SH) 


u gees. 0 0 Cas u 0 
(14 |vl=| 0 ma, 0 eee "0 
w 0 0 mta;; " |lwl |W 


where: @_, are the apparent mass terms, here assumed to be constant 


Assuming no rotation between the fixed earth reference and the system's body 
axis, the ground velocity can be determined. 
Ve =V,+Vy; where Vc, Va, and Vw are the ground velocity, parachute velocity 
relative to the local air mass (airspeed), and wind velocity, respectively. 
(15) V,=V,+V, 
Applying 15 to 14: 


= 


u Mareen 0 0 u 
° . = 5 
(16) Ve=|v}] = 0 [pe (Oin, 0 om v +V,, 
w 0 0 Mm + Ass " lwl [W 
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Equation 16 can be solved numerically to estimate the system response to the 
estimated winds. This estimate can then be compared to the measured velocity of the 
system. Using the measured ground track velocity as the initial wind estimate, the 
modeled ground track was determined. The difference between this modeled ground 
track and the actual measured ground track reflects the errors in the wind estimate. The 
magnitude of these errors indicates the significance of accounting for momentum changes 


caused by changes in the wind. A Simulink® model was implemented (Figure 17). 


Measured Ground 
Track Veloerty 


Vuudot_hat 
dufdt 


Derivative To Wornspacei 


= Va hat MATLAB 
ee Function 


Initial Wind 
Estimate 





Figure 17. Simulink® Realization for Wind Estimation 
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Figure 18 presents the results of this validation. Use of the measured ground 


track velocity as the wind estimate resulted in errors of less than 0.3 feet per second. 
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Figure 18. Wind Estimation Results 


To further assess this technique, the equations for wind estimation utilized with 
hurricane dropsondes® was applied. In this application, the 3-dimensional velocity of a 
dropsonde is determined. The horizontal velocity components of the air mass are then 


found as follows: 
(17) yp @E—A and Vy x P=: 
& & 
To apply equation 17, a Simulink® model was implemented (Figure 19). The 
correction that would be applied to the horizontal velocity is presented in Figure 20. 


These data show that the corrected wind estimate, using the hurricane approach, differs 


from the measured ground track by less than 0.3 feet per second. More errors are seen 
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close to ground impact where additional shears are present. This is not of consequence 


for this study as the flight dynamic data of interest is at much higher altitude. 
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Figure 19. Hurricane Wind Estimation Implementation 
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Figure 20. Hurricane Wind Estimation Results 
These results indicate that the momentum effects can be ignored for wind 


estimation for the selected parachute system. Figure 21 illustrates the comparison of the 
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wind estimate to the winds measured by the RAWIN. Recall that the RAWIN balloon 


was launched only every hour. The closest RAWIN data were used for this comparison. 
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Figure 21. Wind Estimation Compared to RAWIN 


These results demonstrate this technique will provide significantly better 
estimates of winds than using the RAWIN system. By adjusting the weight of the 
calibration system to match the descent rate of the test item, the two parachutes will be 
subjected to the same (as close as possible) atmospheric conditions. Using the measured 
GPS ground track velocities 1s an adequate approximation for wind estimation. Other 
techniques, such as that presented above, may provide some refinement on the wind 
estimates, but the difference is likely to be insignificant for most testing. The key to 
application of this technique is the use of a very stable parachute due to the reductions in 


apparent mass effects resulting from oscillations. 
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DATA REDUCTION 
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Data from all sources were time correlated using the AHRS time tag. As these 


data were recorded at the fastest rate, all other data parameters were correlated with this 
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and density effects were removed from the calibration parachute vertical velocity. 


ficant oscillations 
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The performance of the C-9 parachute is characterized by s 


(Figure 22) in both pitch and roll. Due to the increased mass of the trapped air, the larger 
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G-12 parachute system is not expected to oscillate as severely as this parachute. 
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Figure 22. Parachute Oscillations (video clip) 
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The system instrumentation was all located in the payload and therefore measured 
the movement of the payload induced from the oscillations. The velocity data best 


demonstrates this character (Figure 23). 
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Figure 23. Data Characteristics Due to Oscillations 

To better represent the operational system during this feasibility test, attempts 
were made to modify the C-9 parachute to minimize the oscillations by cutting large 
symmetrical slots in the canopy. This modification proved very stable during tow testing 
as well as the initial test drop. The system was then dropped with one riser extended 
representing a single control input. The objective of this drop was to qualitatively assess 
the amount of drive that could be obtained with this modified parachute. There was 
insufficient indication that this modification allowed enough drive for the guided system. 
Had the results been favorable, additional instrumented drops with the actuator/control 


system would have been accomplished. 


phe, 


Since the trials with the modified parachute did not present the desired results, the 


oscillation in the data needed to be dealt with prior to parameter estimation. These data 


must be corrected to the location of the system's center-of-gravity. 


(18) 


(19) 


(20) 


(21) 


Starting with the position corrections: 
P., =P,+;RL; where P, is the measured position of the payload and L is the 


lever-arm from the location of the instrumentation to the 
center of gravity in body-axis coordinates {B}. 


The velocity correction is then given by: 
gS) 

ee P... = Es oe 

Noting: )R=5RS(@) 


RS cop". Ri (capaiag) 


Substituting (20) into (19): 
aes up 
Vg = P,+3R(@xL)+ i iG 
For rigid bodies, L 1s assumed constant and therefore <L = 0), 
i 


Applying this condition to a parachute did not yield adequate results, as the 


oscillations were still apparent in the velocity data. It is presumed that in an oscillating 


parachute, the forces on the parachute are changing due to the linear acceleration induced 


by the oscillations (apparent mass effects) while the forces on the payload are not 


varying. With this imbalance of forces in the parachute and payload, an effective change 


in center-of-gravity of the system results. 


To illustrate this situation, the lever-arm correction is determined by rearranging 


Equation 21. 


(22) ™ L=ER(Pog ~ P.) 
P, is measured. Pcg 1s estimated using the measured velocity at the payload 
(VL). A low-pass filter is applied to Vy to remove the effect of the oscillations. Using 
this filtered velocity (V_F) and removing the wind velocity (Vw), an estimate of the 
velocity at the center-of-gravity 1s obtained (Vcc). 
(23) Veg *Vip Vy 
This velocity is then integrated from the parachute release point. The estimated 
location of the center-of-gravity was then determined using Equation 22. Figure 24 


presents these results. 
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Figure 24. Movement of Lever Arm Correction 


The mean lever arm corrections are L,=0.0, L,=-0.1, and L,=-3.94 feet. Yavus 
and Cockrell’ demonstrated that acceleration of the air mass and angle of attack of the 


parachute significantly affects the apparent mass coefficients. Attempts were made to 


correlate these data to linear and angular acceleration of the payload. However, no direct 
correlation could be found. Ideally, acceleration of the parachute should be measured in 
future work. With the existing instrumentation suite, it appears the flight test data must 
be filtered to remove the effects of oscillations prior to parameter estimation. Using the 
data obtained from the uncontrolled test drop (311drop1), a filter was derived using the 
MATLAB*® System Identification Tools. A 5“-order low-pass filter with a cut-off 
frequency of 0.009 (in fractions of the Nyquist frequency) was selected. Figure 25 


compares the measured data to the filtered data for velocity. 
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Figure 25. Comparison of Measured and Filtered Velocities 


This filter was applied to the position and velocity data for all drops to estimate 


the motion of the system at the center-of-gravity. The obvious concern about filtering is 


© Registered trademark of Mathworks, Inc, Natick, MA. 


Sy 


the potential of eliminating true system motion while removing "noise." Figure 26 


compares the glide ratio filtered and unfiltered response to a single control input. 
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Figure 26. Assessment of Filtering 


Implementation of the designed filter appears to have the potential of eliminating 
desired response data. The character of the response (magnitude) remains in the filtered 
data but the response time appears to be adversely affected. A filter of this magnitude is 


not desired but may be necessary for parameter estimation. 
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Ke FLIGHT TEST RESULTS 


The control system 1s intended to affect a change in horizontal velocity. This is 
best demonstrated by assessing the glide ratio of the system with the winds removed. 


Figure 27 presents the glide ratio with the measured control inputs. 
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Figure 27. System Response - Glide Ratio 


The results show that a nominal glide ratio of 0.4 to 0.5 exists with no control 
inputs. Potential causes of this induced glide are motion induced by the oscillations, 
imperfections in length of the pressurized actuators, the mathematics of creating a 
horizontal glide ratio which eliminates direction of motion, or errors in the wind estimate. 
This nominal glide ratio does not limit the assessment of the response due to control 


input, as we are interested in the change of glide ratio at the time of control activation. At 
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time zero, all PMAs were pressurized. The system was then allowed to stabilize to a 
‘trim' condition. A change in glide ratio is apparent at approximately 20 seconds with no 
change in the state of the controls. The first incident of change in glide ratio can be 
attributed to the parachute stabilization/inflation process. The remaining data clearly 


show a correlation of glide ratio changes to the activation of the controls. 


A. SINGLE CONTROL INPUT 

Figure 28 isolates the response of a single control input. An increase in glide ratio 
from approximately 0.5 to approximately 1.0 with a time constant of about 4 to 5 seconds 
is observed. The system returns to its oscillatory trim state after about 5 seconds 
following removal of the control input. The reduced magnitude of oscillation or coning 


angle contributes to a reduced rate of descent and increased glide ratio. 
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Figure 28. Single Control Response 
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B. TWO SIMULTANEOUS CONTROL INPUTS 


Recall that the two control inputs can be activated simultaneously. The intent is 
to provide additional resolution (every 45 degrees) in controlling the system. Figure 29 


presents the response with two simultaneous inputs. 
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Figure 29. Two Control Response 


As exemplified by this figure, there is ee in performance with two 
control inputs over that achieved with one. In fact, the data indicates reduced response 
results from the simultaneous activation of two PMAs. This reduced performance 1s 
likely due to leading edge collapse (as observed in the ground-to-air video, Figure 30) of 
the parachute with two control inputs. The magnitude of the oscillations is not reduced as 


dramatically as with a single control input. 
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Figure 30. Parachute Collapse from Two Control Inputs (video clip) 


The vertical velocity of the system was observed to decrease with a single control 


input but remained essentially constant with two simultaneous control inputs. This 
reduction is likely due to the reduced "coning" angle from the oscillation damping effect 


of the control system. Figure 31 presents the vertical velocity with control actuation from 


a single and a double input. 
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Figure 31. Vertical Velocity Response to Control Input 


C. ATTITUDE RESPONSE 


For control system design, understanding the heading and yaw rate responses of 
the system is critical. The heading and attitude rates, as measured by the AHRS, were 
correlated to the control inputs for analysis. Actuation of the control system resulted in 
significantly decreased parachute oscillations. Figure 32 illustrates the roll rate data 
collected from the AHRS. Whenever a control was activated, the attitude rates were 


significantly reduced. 
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Figure 32. Roll Rate Response 


The root mean square (RMS) of the roll rate reduced from 16.8 degrees per 
second down to 2.1 degrees per second with control activation. This behavior was 


observed each time a control was activated. 


A plot of the total displacement (or coning) angle from vertical also demonstrates 


the reduction in oscillations. 
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Figure 33. Oscillation Angle with Control Input 
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Figure 34 presents heading, measured by AHRS, and the control activation. 


Figure 34. Heading Response 
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With careful examination of these data, one can correlate changes in heading with 
control inputs as indicated with the arrows on the graph. However, the response appears 
to vary greatly. The magnitude and direction of heading changes vary with the activation 
of the same control. The cause of this variation in heading response has not yet been 
identified. Factors such as the coning angle at the time of control input must be 
evaluated. | 

During initial modeling efforts, it was apparent that the rotation rate of the C-9 
parachute severely limited the efficiency of this control concept. Excessive control 
inputs would be needed to affect control in a single cardinal direction. Therefore, it was 
crucial to this study to determine if the operational G-12 system exhibited similar 
characteristics. The results from this demonstration clearly show a significantly reduced 


rotation in the larger G-12 parachute as compared to the C-9 parachute. The C-9 


rotations can be observed in Figure 35. 
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Figure 35. C-9 Rotation Rate (video clip) 
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igure 36 presents a video from a camera installed on the payload. With the sun 
as a reference, the very low rotation rate of this system can be observed. A ground based 


video (Figure 37) also demonstrates a low rotation rate for the G-12 parachute. 








Figure 37. G-12 Rotation Rate (video clip) 
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VI. NAVIGATION SENSOR MODELS 


The AGAS is expected to include two navigation sensors: 1) a commercial Global 
Positioning System (GPS) receiver for position determination and 2) a heading reference 
assumed to be a magnetic compass for this study. To assess the effects of navigation 


sensor errors, models for each sensor must be incorporated into the simulation. 


A. GLOBAL POSITIONING SYSTEM (GPS) 


GPS error sources include errors induced by the atmosphere (ionospheric and 
tropospheric), multi-path, receiver noise, satellite clock noise, and Selective Availability. 
Modeling techniques for GPS range errors resulting from these sources have been 
developed and validated® to model range errors and not errors in a Cartesian reference as 
desired here. Cartesian (x, y, z) errors would therefore have to be formed from the range 
errors for implementation in this simulation. This necessitates application of a numerical 
solution like maximum likelihood techniques.” Although this implementation is 
relatively trivial, the computation resources required severely limit the simulation speed 


on a Personal Computer (PC). Therefore, a variation of this approach was implemented. 


1. Selective Availability 


Selective Availability is a means of intentionally inducing errors into the GPS 
satellite signal. The DoD induces these errors to restrict use of the full precision of GPS 
to unauthorized users. Authorized users must apply a receiver capable of processing the 
cryptographic codes to remove these induced errors. Although the AGAS concept could 


incorporate an authorized receiver, it is desired to utilize a commercial GPS receiver for 


cost savings. Also, with airdrop the loads may not be fully recoverable and loss of the 
authorized receivers would not be desirable. 

The errors resulting from Selective Availability are not stochastic in nature. 
Therefore system identification methods were employed to obtain a reasonable error 
model. Data were collected at the Yuma Proving Ground Satellite Reference Station. An 
unauthorized GPS receiver was placed on a known survey point. GPS position data were 
collected for over two hours. The position data in the three Cartesian axes were 
differenced with the surveyed coordinates resulting in the Cartesian errors. These errors 
represent all GPS error sources identified above. Figure 38 illustrates the apparent 


random nature of these data. 
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Figure 38. Measured Unauthorized GPS Errors 


To obtain a model of these data, the MATLAB* system identification toolbox 


was utilized. An ARMAX"® model was utilized with the input being white noise and the 
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output being the position errors shown above. The ARMAX model incorporates a 
prediction error method with a model represented by a set of difference equations of the 
form: 


(24) A(q) y(t) = B(q)u(t — nk) + C(q)e(t) ; where y and u are the outputs and inputs of 
the system, respectively. 


The coefficients A, B, and C are polynomials that describe the model's difference 
equations. The prediction error is minimized using an iterative Gauss-Newton algorithm. 
The ARMAX function returns a matrix of the polynomial coefficients. This matrix, 
referred to as THETA format, can then be transformed into a transfer function using the 
MATLAB® command TH2TF. This technique resulted in the following transfer function 


that was utilized in the overall system model to obtain GPS errors: 


z* -1.5302z* + 0.2608z* + 0.2566z+ 0.0192 


(23) 00 
z -2.6500z° +1.9582z° + 0.0337z -0.3420 
Figure 39 presents the output of the GPS error model including Selective 
Availability Errors. 
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Figure 39. Modeled GPS (Unauthorized) Errors 
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The transfer function input 1s white noise initiated with a random seed ensuring 
variable errors are introduced from simulation to simulation. To assess the adequacy of 
this model, the mean, standard deviation, and root mean square were calculated for the 
measured and modeled GPS errors. The sample of measured errors presented above has 
a mean value of approximately zero meters in each axis and a standard deviation of 17.2, 
28.8, and 21.1 meters in the x-, y-, and z-axes, respectively. The modeled results 
demonstrated mean errors of three to six meters with standard deviations ranging from 25 
to 35 meters. The root mean square errors for the model were found to range from 26 to 
37 meters for three independent simulations. The model produces a reasonable 


representation of the measured GPS data. 


dy Model Without Selective Availability Errors 


With Selective Availability turned off, that is no induced errors, a commercial 
GPS receiver is capable of navigating with greater accuracy. A GPS Error Model was 
derived considering a noise structure proposed by Draper Laboratory.'’ This report 
models a P-code GPS receiver incorporated into the Honeywell Embedded GPS/Inertial 
Navigation System. The noise model incorporates two components, accuracy and jitter. 
The accuracy noise component is considered exponentially correlated noise. The jitter 
component consists of two elements: an exponentially correlated noise component with a 
faster time constant than the accuracy component and a uniform uncorrelated noise 
component. The GPS position noise model, suggested for GPS-only operations (no 


inertial aiding), was adapted for a commercial grade (C/A-code) receiver by adjusting the 
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accuracy and jitter standard deviation specification. Table 3 presents the original and 


adapted models. 






Draper P-code C/A-code 
30 45 

Tier ae a ae | ——¢ —— 
WO sar@) | T*sqnG) 


Note: q is the quantization interval 
Table 3. GPS Position Error Model 














This model was incorporated into the Simulink® simulation. Figures 40 and 41 
illustrate the results obtained from this model for a 400-second simulation and i acl 
simulation, respectively. The standard deviation of the 3-axis error for this simulation 
was 56.3 feet (17.2 meters) which is close to the specifications for an commercial 


receiver with Selective Availability off. 
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Figure 40. GPS Error Model - SA Off 
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eer 41. GPS Error r Model - SA Off f (50- coal | Simulation) 


B. HEADING SENSOR 


The heading sensor is assumed to be a magnetic compass for this study. Two 
components of errors are considered here, a static error or bias and a dynamic (noise) 
component. System specifications for the Attitude Heading Reference System (AHRS) 
provide a static error of +2 degrees (+ 1 degree with velocity aiding) and a dynamic 
component of +2 percent. The AHRS incorporates rate-gyros to obtain 3-axis attitude 
rates and attitude data. Specification sheets of a low-cost digital magnetic compass 
produced by KVH Industries presented similar accuracy statements. The static error 1s 
incorporated as a bias element in the Simulink® model and is set as a uniform random 
variable at the start of each simulation. The dynamic component is found by adding 2% 
of the current heading reading. Figure 42 presents a 400-second simulation of the 


heading error. 
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Figure 42. Modeled Heading Error 
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VU. EQUATIONS OF MOTION 


The kinematic equations of motion are determined assuming the parachute 
and payload system are a rigid body. A six-degree of freedom model is considered with 
the motion about the center of gravity. The notation and derivation follow the convention 
presented by Isaac Kaminer’’ as well 5 by Antonio Pascoal and Carlos Silvestre’. The 
derivation first considers an origin displaced from the system's center of gravity by a 
distance P,. In the implementation of these equations to this problem, the origin is 


moved to the center of gravity of the system. 


A. NOTATION 


The following illustrates the reference system used for this effort: 





Figure 43. Coordinate Convention 


The following notations are used in the development of the equations of 


motion. 


Capital letters/symbols denote vectors or matrices 

Small letters/symbols denote scalars 

{B}: name of the coordinate system, B denotes body-axis, W denotes wind- 
axis, and U denotes inertial axis. 


- ©R: rotation matrix from coordinate system {4} to coordinate system {C} 

- °P.: position of point D, measured in {C} and expressed in {C} 

- “(°P,): position of point D, measured in {C} and expressed in {4} where 
“(CPp) =é RCPp Py 

- °V,: velocity of point D, measured in {C} and expressed in {C} 

- ©Q,: angular velocity of point D, measured in {C} and expressed in {C} 


2 . : time derivatives in the body-axis {B} 


- ('): time derivatives in the inertial-axis {U} 


NOTE: if the symbol for the coordinate axis 1s omitted, the inertial-axis 
{U} is assumed. 


B. ASSUMPTIONS 


Throughout this effort, the following assumptions were considered: 
- Rigid-body system 
- Only the motion after complete parachute deployment is considered. 


- Non-rotating earth, 1.e. ECEF coordinate equals inertial coordinate 
- Wind axis equals inertial axis 


C. DERIVATION OF EQUATIONS OF MOTION 
Starting with Newton's Law, F = ma 
Vike 
The inertial velocity of the body is: 
“V =" RP V, + QxP, S ; where Po is the position of the parachute 


reference to the origin (center-of-gravity). 


Sy 


Taking the derivative: 
, d ; 
Come ie RO, +QxP,+QxP, {eR('Y, +Qx P,); where P, =0 due to 
rigid body assumption 
recall: 


Pp 
Gime ~R=-R S@O.), *O,=\¢ 


QO =r @g@ 
and S’?Q,)=| r 0 -—p| (skewsymmetric matrix) 
-q p 0 


Substituting (27) into (26): 


(28) Py =$R)SAV, + xP b+GR-SPM4)(V, +2%Po) 


Noting: 
(29) SGQ,) Ve OF ae. 


Substituting (29) into (28): 
Uy? U d B : UnpnsB B 
rs -$R\ 4 OP, |+SR Q,x(2V, +QxP,) 


Rearranging Terms: 


Ve SRY, - P, + 5R Ox eax 2+7Y, | 


Consider forces in the body axis: 
®F=f$R°F=$R-M,"V,; where: M, =ml, 


iss 


= hae OlesR Q,x{- Py x QV, 


=M_ER|5R| 


aM ov, —P, Oh +, > Q x(- je x 2477, | 


(30) = M,<'V, -M,P, x Q47.0,x(-M,P, xQ4+M, aye) 


Again applying the notation used by Pascoal by defining the general mass matrix: 


31) Vix Me Me 7 il —m’ P, x 
M iv m”P, x JG 


Equation (30) reduces to: 
d B = B B 
F,=M,—'V; + Mp Q+? OQ, x (Mig Q+Mz ?V,) 


Solving for velocity in the body axis: 


(32) “*Y, = Vie Ono x Mi Q +M,°V; Ni M, Fs 
Now applying Euler's Law: 
(33) “%L,="N, -M,P,x"V,; angular momentum equals total applied moment; 
where: Af = ml,; [3 is a 3x3 identity matrix and m 1s the system's mass 
Recall: 
VARY, 


LY 3 GRAV ERY, 
v RPV, =" R(Ox®V, } 
Oey U d B U B | 
34 °V,=5Ro V,+UR(Qx?Y, ) 
Substituting (34) into (33): 


"La='N, -M,P, x & RE*V, +ER(Ox! V,) 
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(35) 


(36) 


(37) 


(38) 


d 


oh ay ad ie cf Su 
fey es |: es dt 


LGR L=GRO*L+ SRL) 


at 
di 


y REAL+SROMLEN, -M,P, “( UR 


®V7, +" RIOx®V,} 
By definition: 
FL=l 32; where Ip is the system's moment of inertia matrix 
Hie HAG 
Substituting (36) into (35) 


URI, Q42R(Qx1,Q)="N, -M,P, x v Ro, +2R(Ox? V,) 


Rearranging Terms: 
uN =" RT,04+9R(Qx 1,0) +MpP, “( Z REV, +R V, ) 


Transforming into the body axis {B}: 


2N = PRN, =1, <2 +(Qx1,Q)+M,P, (or, +(Ox?V, ) 


*N,=1, “0+ (Qx 1,Q)+M,P, x(Qx?V,)+M;,P, x SV, | 
Applying (31) to (37): 


BN, =M,20+Mq{ 7; }+OxM,0 + Mar (Ox? V, } 


Rearranging: 


d 


Ogee : 
Gas Mar $27, + Oxy + My (OXY) +M, os 


dt 


Combining (38) and (32) in state space form: 


5) 


da : 
| 3 | M, —°V, + Mg Q+?.Q, x IM ,2+M, °V,} 


NG 


d 
M, 52+ Mre| SV, |+OxMO+Mg (Ax? VY, | 


d ; . 
<V, M, M2 Q+?Q, x {M2 2+M, °V, i Meas Oe aa 

(39) ST d dite 
Cece Man( G's |+ 2x Mp My (0x? Y,) 0 Bie 
at 


D. APPARENT MASS 


As a body accelerates through a fluid, the fluid itself must accelerate to 
accommodate the motion of the body. Resultant forces and moments are applied to the 
body. A common method for accounting for these forces is to include added or apparent 
mass terms in the equations of motion. Sir Horace Lamb performed the onginal work on 
the effects of accelerating fluid on a body."* Lamb derives the apparent mass effects and 
identifies 15 independent terms. Cockrell and Doherr’” showed that for a revolution of 
the body about a plane of symmetry, the number of independent terms reduces to the 


following: 


@,, - motion along the x-axis 
a,, - motion along the y-axis 
a, - motion along the z-axis 
a,, - motion about the x-axis 
a, - motion about the y-axis 
a,, - coupling between motion along the x-axis and about the y-axis. 
a,, - Coupling between motion along the y-axis and about the x-axis 


These terms have been demonstrated to change significantly with 


acceleration of the body.'° However, estimation of these changes is very difficult and 
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beyond the scope of this research. Therefore, only the scalar values of apparent mass and 


moment of inertia terms are considered. 


I Estimation of Apparent Mass Terms 


Doherr and Saliaris'’ demonstrated that for an equivalent arachute, four 
q Pp 


independent apparent mass terms exist: 


(40) A, = Ga(H*)%} ; D, is the profile diameter of the parachute 
3 


1] 
(41) G); = A222 ~ 5 33 


(42) Cr -O—) 2Oe 


Po) ee —c.. = 0.0480), a, 
Z, Effects on the Equations of Motion 


The force imparted on the body by the air mass can be determined by starting 
with the equations of motion for the body itself (Equation 39) and replacing the mass 


matrices with apparent mass matrices as shown here: 


Mr D Ar Mrr D> Arr: Mar > -ArrR: Mr D Ap 


: M,<V, + Mp Q+? OQ, XM pQ+M, *V, } 
(44) f |- . ? 
N; My 2O+Myy[ ZV, | + 2xM OM gy (O%'Y,) 
i ae Ga 
(43) A-=| 0 a, O01}; A,=| 0 a, 
0 0 @, 0 0 
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o> OO O Gane OS One 
(46) Arp =| Go, ~O-RONRA., —|C,. MOR ae 0 O 
O07 OG SE: OC) 0 Or 


(47) Ae = = Ae. 


RT 


Substituting (46 and 45) into (44): 


ad ; 
77 Ay —7V 5 + Ag+? Qy x 1A,,2+ A, V5} 
Cay BN m d d 

y Ay ZO+ Ap| SAVy ) +x A,Q+ Apr (2x? V,) 


In addition to those forces and moments presented here, an additional term is 


required in the moment equation to account for the additional moments generated by 


changes in momentum due to air acting on the parachute. The term *V x A,°V is added 
to (48). 


d 
a | A, 'V; =p AON Ox {Ang Q + A, a, 


(49) i : 
°N Ay T+ Age SAV | +2 Ag+ g(t Vy AV x dV 


Combining (49) and (39) and rearranging terms: 


< ay, Mi M1, QQ, x {Mi,Q+M'2y, tt tg Ler. 
= M, 
49 usta £20 , |+ xm, Q4+M.,( ONE i VxA, ay |") 9 Mi Bal 
dt 
where 


va Mi) {Me Mo * dn 
Vee eee eA. 
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E. EXTERNAL FORCES AND MOMENTS 
The external forces and moments on the body include the effects of aerodynamics 
(including actuation of the control surfaces) and gravity. 
| Eee pe ee 
?N, = Nero t Nay + Noriv 


18 C . 
|” are not included as the convention used 


Buoyancy terms as discussed by Gocke 
herein does not include the mass of the trapped air in the overall systems mass. As 


Gockel demonstrates, if one is to include the mass of the trapped air a compensating 


"buoyancy" term must be added to the equations of motion. 


1. Aerodynamic Forces and Moments 


Using a first order Taylor series expansion about a nominal trim condition and 


following Kaminer's’’ notation, the aerodynamic effects are found to be: 
Popp = + Or Or.x FOR 
N zero = N, + ON, +ON,X+O0N,A 
where x and x are the state vector and first derivative of the state vector and Ais 
the control vector. OF and 6M are the aia derivatives of the forces and moments with 


respect to each parameter in the state vector and their first time derivative. 


Non-dimensionalizing and putting in state space form: 
F ysro = WSo(Cro + OC pr, tC 5X +5C pA} 
N sero = 95,4, (Cy, tOCy, tOCy x +0Cy, A} 


where: 


=o 


Cr = (oa Cy Cee ee ie Gel 


Cr, 7 ies Oo CF oe oy Call 
Cy Cy Cy, Cy Cy es 
6C re =||Cru Cy Cy Cy, Cy, Cy, 
Cx, Cy Cay Co, Cy, Cz 
oF Oss oF ee L@ ee 
OC =| CG er, CF Cu Cur 
Cw Cz Cw Ce Cu Cy 
C xs Cy Cw Cy Cy Cy 
OC rm =|| Cy Cy Cry Cy, Cy, Cr 
Cy Ca Cy Cy Cy Co 
oF Crs Ch Ci C1; Cy 
O Cvs =}! Cary Cie Cp Cup Cus Cu 
Cn Cn er Cx Cng Cn 


OC p= len Cw Gr OC yy = (om Cun Cul) 


De Gravity 


The gravitational force acting on the body in the inertial reference frame is given 
by: 
0 
oy 
Fopsy =| 0 


mg 
The force in the body axis is found using an Euler angle transformation: 


Q 
ey ye na 0 
mg 


3 
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cosy cos@ siny cos@ —sind 


where ~R=|cosysin@sing-sinycos¢ sinysin@sing+cosycos¢? cos@sing 


cosysin# cos¢+sinysin@ sinysin@cos¢—cosysing cos@cos¢ 


Recall that P, is defined as the distance from the origin to the center of gravity. 


The gravitational force creates a moment about the defined origin. 


0 
Nee = Po hee, = eh OM iatemes willbe dened: 
mg 
3; Total External Forces and Moments 


The total external forces and moments are now given as: 


0 
PE = qS,\8Cr, t+ 0Cp, tC p,X+OCp, AMER) 0 
mg 
m0 
Nz =]S,d, Cy, +OCy, OC y .X+OCy AS+P,FR\ 0 
mg 


In matrix form: 


0 
JemlgS f.00 0 Cr, FOC; FOC, soe, A, Ae] 
®N,| | 0 GS,d4,f; | Cy, tOCy, +OCy,X+5Cy AS] | PLR 


oO 


mgs 


| Ve COMPLETE EQUATIONS OF MOTION 


This derivation assumes that the system is a single rigid body and resulted in a 
general set of equations of motion. These equations could then be applied to the 


parachute and payload separately. For example, to fully capture the oscillations of the 
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parachute system, it may be beneficial to establish a coordinate system at the parachute 
aerodynamic center. The model would then incorporate the applicable components of 
these equations for both the parachute and payload. The general equations of motion are 


completely defined as: 


a My Mig Q#? Og x Mig 2+ Mr "V5 j 
ee — My" {Min( LAV, }+ Ox MpQ+ Miy(Oo°Y,) 4° x4 27 eene 
dt 


0 


0 Mil 0  GS,d,1, | (Cy, +oCy, +OCy,%+5Cy, AS] | PER 


Oo 


ie 0 ie 0 hema 
where: 


wi Mn) | Me*4 Ma * 
M pr M, Mls Sp ile M , Ap 
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Vill. ANALYSIS OF EQUATIONS OF MOTIONS 


A. EQUATIONS OF MOTION EXPANSION 


Expansion of the equations of motion is necessary for parameter estimation. This 
expansion also serves to validate the derivation of these equations against existing works 
in the field. First, the coordinate system axis 1s selected aligned with the parachute body 
axis but displaced from the center of gravity by a distance z. Expanding the equations of 


motion: 


xX =(n+a@,, )u+(mz+a,,)¢+(m—a,,)pr—(m+a., rv +("+a,,)qw 
Y=(M+Qy)V—(mz—-Q,,)p+(m+aQ, )ru + (mz +@;)gr —(m+a;;)pw 

Z =(m+@3,)W—(mz +@5)q° — (m+ a,,)qu—(m-a,,)p° + (M+ ay) pV 

L=(1, + Q4,)p — (mz -a@,,v+ UZ, — 1, — Ass )gr — (mz + G5 )ru — (mz — A2,) pw + (Q,, — As; UW 
M =(1, +@55)9+(mz+a,5)u+U, -—1,+@,,)pr — (mz +a,;)(qw-1v)- (a, —@;3 uw 

N= 


These equations are validated by the work of Cockrell and Doherr””. 


Taking terms like rv, qw, etc. to be small and eliminating non-linear terms like p*: 


Ee — (ae ze, 
Y =(7+a,, )v—(nz—a,,)p 
Z =(m+aQ,,)W 
L=(,+@,,)p—(mz-Q,,)V 
M = (1, +@5)q + (mz + @,,)u 
N=I sr 


Forming forces and moments as an initial (trim) value plus a perturbation and 


applying gravitational forces and moments. 


X = AX —mgsin(@) =(m+a@,, )u+(mz+@,5)q 
Y = AY —mgcos(@)sin(¢) = (m+a@,, )v—(mz-a,,)p 


Z=Z,+AZ+megcos(@)cos(¢) =(m+a,,)w 

L = AL + mez sin(@) = U7, +@4,)p — (mz —a,,)v 
M = AM - mgzsin(@) = 1, +a@55)q+(mz+a,,)u 
N=AN=I1,r 


Subtracting initial conditions as defined as Xo=0, Zo=-mg, YO=0, Lo=0, Mo=0, 
No=0, 89=0, ®o=0 and applying small angle assumption (e.g. cos@ = 1,sin @ ~ @ ) two sets 


(Longitudinal and Lateral-Directional) of uncoupled equations are formed. 


B. LONGITUDINAL EQUATIONS 


(m+a,,)u+Gnz+a,,)¢g =AX —mge 
(m+a,,)w= AZ 
(mz+a@,,)u+(1, +a5,)q = AM —mgz@ 


The external forces and moments are approximated with a Taylor-Series © 


expansion: 


W 


0 0 0 


AX = Akg oe eee le eg Negras 
W. W, q q 


AV AZ. AZ, = + AZ a AZ. 4 AY Gg EA Gee 
W W. W. 3 a 


0 0 


AM = AM, -——+AM,——+AM, —-+AM, —-+AM.q+AM,g+AM,6 
W W. W. W. 


0 


Reducing terms that have no apparent effect on forces or moments in a particular 


direction (e.g. AX, w, AZ.) and substituting: 


"+ AX, ——mgO+AX,d 


u 
LUO eI Resear ES ler Le) 120% WV, W, 


0 


W Ww Ww 
om + ass) oF = AAW ay, Oe eae 


0 
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W, (mz + a5) — pats +Q..)q = AM, 7M, oS 7+! wt M9 mgz@ + AM,6 


Defining the stability derivatives: 


a AX, | = AZ, AM, 
W (m+a,,) W,(m+a53) WU, + ss) 


where: i=u,u,w,w 


_ es ae L—— = j 
, (m+ay)- (arcs) ’ (1, +55) 


where: j/=4,g,0 
Establishing the coordinate system now at the center of mass of the system, 


substituting the stability derivatives and rearranging terms: 








(‘- v ae Qs g=X,—- mgé@ 4, 450 
W, Pelee) W, W(mt+a,) W, 
Z 
(| Sa ee 
W, sal Wy OW, 


a u 
W,| ———- -M, |— -W.M, ea, = W. madd iM +M.6 


In state space form with 6 =q: 











a ut —mg u xX 

1-X, 0 2 C= a 0 0 ee 5 

( ) W(m+a,,) W, : W,(m+a,,) W, W, 

0 (1-Z,,) 0 w Zs wl |Z, 
—/|-= 0 Zz. — 0 —- i rs) 

Qis _ -WM l 0 Wy W, W, W, 

(+55) ae q| |WM, WM, M, 0 q| 1M, 

0 0 1|L @ 0 0 1 0 G 0 


c; LATERAL-DIRECTIONAL EQUATIONS 
(m+ hu — (mz - a,,)p = AY —mg¢ 


—(mz—-a,,)v+(/,+a,,)p =AL+mez¢ 
Il r=AN 
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The external forces and moments are approximated with a Taylor-Series 


expansion: 


AY = AY, + AY, +AY,p+AY,p+A¥,r + AY? + AY,6 


0 0 


AL = AL, —— + AL,—+ AL, p+AL,p+AL,r + AL,? + AL,6 
Wi We ‘ 


AN = AN, ——+AN,——+AN, p+AN,p+AN,r+AN,F+AN,6 
YW "W Pp Pp r r 


0 0 


Reducing terms that have no apparent effect on forces or moments in a particular 


direction (e.g. AY, p; AN, ) and substituting: 


Vv V v 
W (m+a..)—-—(mz-—a,,)p = AY. —+ AY, —+AYr+AY.6 —m 
S(O Gas) Ge (ts Cs) P= AY, G+ AY, + AY 7 + AY, — mgp 


Oo 


v . v v 
Woke — Ora) a + Ue Ry Sik” eee a OD 


0 


fF =AN, —+AN ,p+AN,r+AN,6 
Ww. 


Defining the stability derivatives: 


AY, AL, AN. 


z z 


— 


Y=; L, = _; N= 
W (m+a.) Wd +o) Wl 


al 


where: j=v,v 


SS L. =~ } Ne - 
(m+ G2.) (1, +@4) 


i 


Ne AL, AN, 
ig 


where: j= 7p,p,r,7,6 © 
Establishing the coordinate system now at the center of mass of the system, 


substituting the stability derivatives and rearranging terms: 
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(= 2 ae ee es _ales a 
W, W. (mm 5 +n) W, Wo W W,(m ae oo) 
Om v 
| a = W oL,— +L p+lr+l.o 
oF +Qy, 5) p= W, pP r } 
PAWN PN pt N rAN 
0 
In state space form with ¢ = p: 
(1-¥,) i ee 0 0 v yi 0 . = mg a Y, 
W,(m+a) WV, v W, Wy (m + a2) W, W, 
ip ( Sy 1 0 0) QA ee 0 p \+|L; 6 
se ; rol #1 [eM M,N, 0 r | |N, 
4 0 0 11 d 0 1 0 0 d 0 
D. SUMMARY OF EQUATIONS 


In summary, two sets of independent equations have been derived for the 


longitudinal and the lateral-directional motion. These equations form the models to be 


used for parameter estimation. 


Longitudinal: 
ite 0 
0 1-Z 
= u a WM, 
(1, +@s5) 
0 0 


Lateral-Directional: 


(3a) a2 
: W(m+ax 
a Capes a l 
eet ns 
0 0 
0 0 























Qs a XY 0 x, —mg ue X 5 
W,(m+a,, W, : W. Wi(mtay,)|M%\ |, 
0 we Bian - a ; w Fe a 5 
1 0 W 3 W. W, WV, 
q WM, WM, M, 0 q M; 
0 1iL @ 0 0 l 0 g 0 
of %] 7 y te % ___mg NOY aac 
W, " W, W, W(mta.) |] % We 
OO 2) | La ee 0 fe) Ae ees 
Pyle WANES AN NS 0 r N, 
Ogi rio 0 d 0 
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Non-dimensionalizing: 
Ff =q8S,C, and M=GS,d,C,,; where: q, So, and do are the dynamic pressure, 


parachute reference area, and the parachute reference diameter, respectively. 



























































Longitudinal: 
qs Qs 
l- ee 0 0 
(m+ a1) is W.(m+a,) = 
0 poe eg 0 0 if 
(eee Sila 
W Dis > qo, 4, Cap = W,q5,d, @ . ] 0 4 
°| UL, +@s5) dT, Hom) I 6 
| 0 0 0 
ite —m 
qs, 0 0 0 Con 0 0 — 2 u Ce 
(m+a,,) C 75, Ww. W 
7 0 ° 
| ee 0 oll} 0 ag cere Py, || hice 
= (m+a3;) 2W,) W, ral lee 
= 0 ° 
0 0 Bodo _ 0 W Cig cele ah Cc q Cus 
as +.) 0 u 0 w 2W, Mq 9 0 
L QO. 0 0 ] | O O iI @) - =i all 
Lateral-Directional: 
(eee 2 olf v 
(m+a,.) q5W, WwW 
a,,-Gs.a C,. - 
W( 24 q a 6 =v) ] 0 0 P = 
1. +a, c 
0 0 1 0 j 
| 0 0 01 
GS, 0 0 0 C d, Cy, d, Cy, _ Ng 
(m + a5, “ 2W,) W, 2W,) Wy qsW, |) _¥_ Cys 
G54. d d W,| | W, 
0 0 OHI WC (0 ae 0 
= (I, +a,,) eas ca os 4) =e P \+| Cis |6 
0 0 Bt Olle, (4c omic PAP 
Ls, | o~Ny 2W, Np 2W,, Nr g 0 
Y HE l 0 0 


Which are in the form /,x= 4.x+Bu; define: A=/,'A, and B=1,'B, to get 


in the form: x = Ax + Bu 
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IX. PARAMETER ESTIMATION 


A. OVERVIEW 


Parameter estimation is the process of determining aerodynamic stability and 
control derivatives from flight test data. A large number of techniques have been 
employed for parameter estimation in aircraft research. This study investigated the 
application of Maximum Likelihood Estimation methods described by Richard Maine 
and Kenneth Iliff of NASA”’. 

Wind tunnel data on the stability and control derivatives for this system are not 
available. This effort encompassed an initial start at tackling this problem. Parachute 
aerodynamics are truly non-linear and significantly coupled across axes. The maximum 


likelihood techniques implement linear models with no cross-coupling. 


B. MODIFIED MAXIMUM LIKELIHOOD ESTIMATION 


Once the system is modeled, the system's equations of motion, expressed in state 
space form, are utilized to model the dynamic response due to a simulated input with an 
initial "guess" of the aerodynamic parameters. The modeled output is compared to the 
flight test results; the unknown parameters are adjusted as to minimize the output error. 
Detailed descriptions of the MMLE3 algorithm as implemented in the MATLAB are 
provided in the user's guide”. A block diagram of this concept is provided in Figure 44. 

A coning motion characterizes the dynamics of a flat circular parachute. This 


motion will likely not be identified during the estimation process due to the averaging 


© MATLAB is a registered trademark of Mathworks, Inc., Natick, MA 
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techniques implemented. Karl Doherr suggests” the addition of noise to the average 


normal force or pitching moment coefficients to compensate for this. 


Mathematical 
Model 
Maximum Likelihood 
Estimation 


Figure 44. Parameter Estimation Concept 





Ge RESULTS 


The MMLES3 algorithm as implemented in MATLAB386 (an obsolete version of 
MATLAB*®) was successfully run for several test cases. First, MATLAB386 was hosted 
on a Windows 95° on a PC. The program needed to be mun in the DOS mode and could 
not be run in a Windows shell. In addition, several memory management conflicts had to 
be resolved before being able to run the program. After successfully launching the 
program, the Naval Postgraduate School's (NPS) MMLE3 implementation for parameter 
estimation was executed. Several test cases were executed including the T37 and F14 
aircraft simulations. In all cases, duplicate results to those obtained by Graham** were 


obtained. 


© Registered Trademark of the Microsoft Corporation 
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The NPS implementation was then modified to incorporate the equations of 
motion for the parachute (identified above). Simulated data (using estimates for the 
stability derivatives) were generated and qualitatively observed to have the same general 
character of the flight test data. However when utilizing the simulated data to estimate 
the applicable stability derivatives, the algorithms failed to predict the stability 
derivatives utilized to generate the simulated data. The program's execution would halt 
after receiving errors indicating the HES matrix was near singular. This matrix rotates 
and scales the gradient to provide a single step convergence for the quadratic cost 
function.*” Initially, it was believed that this error resulted from the model being 
overparameterized. After significant investigation and numerous trials using simulated 
and flight-test data, no results could be obtained using these algorithms. 

The implementation of the F-14 models was then analyzed again. The F-14 
equations of motion were utilized and the F-14 parameters were modified to reflect a 
slower aircraft (approximately the speed of the descent rate of the parachute). Simulated 
data were successfully generated but when applying the MMLES3 algorithm, the HES 
matrix again went singular. It was concluded that the MMLES3 algorithm, as 
implemented in MATLAB386, could not provide estimates of parachute stability 
derivatives. There appears to be a numerical sensitivity in the algorithms to the relatively 
low numbers representing the parachute system. Further work may result in the 
successful application of this algorithm to the parachute problem. However, this 1s 
beyond the scope of this study. It is also recognized that the utility of results, should they 
be achieved, may not be sufficient. This implementation requires linearization of the 


equations of motion. The parachute dynamic behavior is clearly not linear. In addition, 


other research”® has demonstrated that the stability derivatives themselves are functions 


of angle of attack, and may even vary with time. Figure 45 shows how the drag, lift, and 


moment coefficients change for a typical flat-circular parachute. 
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X. SYSTEM MODEL 
Without successful MMLE results and in the absence of sufficient wind tunnel 
data for the C-9 parachute, a point-mass system was assumed with the only forces on the 
system being drag and weight. The applicable equations of motion are: 
(50) X =(m+@,,)u =—Dcosy cosy 
uy) Y=(m+a,,)v=—Deosy siny 
(52) Z=(m+a,,)w=-Dsiny+W 
where: W is the calibration system weight, y and y are the flight wath angle and 
yaw angle respectfully. D=qC,S D is drag. 


The reference angles are 


(53) ° Ww lies —-WwW , u 
re SS eS OO ie TT 
Ve V IV, ee w- 7 - w~* 
Substituting, rearranging terms, and putting in state space form: 


=| 


u AO 0 0 Cos u 0 
(54) S|) ee M+ Qo 0 Sa Daa | I 
Ww 0 0 M+; 4 wl |W 


(2) where: a, are the apparent mass terms, here es to be constant. 
Noting that: 
Ve =V,+V,; where Vc = ground velocity, Va=velocity relative to the local air 
mass (airspeed), Vw=wind velocity. Assumes no rotation between the fixed earth 


reference and the system's body axis. 


ey 


Equation (54) presents the equation for V,. Substituting, the EOM for estimating 


ground speed results. 


=] 


u tg alo ae 0 0 Cs 
(56) ~~ V..=| ee ee Te 0 a 0 |++V, 
wi. 0 0 m+a;, " lwl IW 


Equation (56) can be solved numerically to estimate the system response. To 
estimate heading, a constant rate of rotation was assumed. The flight test results for the 
C-9 parachute showed significant rotations and changes in rotation rate with control 
activation. Flight testing of the larger G-12 parachute showed a mean rotation rate of 
approximately 1.8 degrees per second. A standard deviation of 1.0 degree per second 
was assumed based on qualitative observations during flight test. An insufficient sample 
size was available to accurately determine a standard deviation for rotation rate. A 
normal random number generator at the start of each simulation determines the rotation 
rate. That rate is then integrated to provide heading. This model was incorporated into 
Simulink®. Multiple wind profiles were incorporated into the simulation using a lookup 
table based on the current system altitude. The equations of motion are contained in a 
MATLAB*® function block which Sais script file titled DOF3. This script file 
(Appendix E) also incorporates the parachute system parameters and apparent mass 
terms. 

Figure 47 presents the measured velocity data from flight test as compared to the 
modeled velocity data for an uncontrolled drop. The 'noise' in the measured data results 
from the velocity being measured at the payload which 1s experiencing significant 


oscillations. Since the point mass model does not incorporate these oscillations, no 
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‘noise’ is apparent in the modeled data. As is demonstrated in the graph, the velocity data 
agree very well for this uncontrolled condition. For this run, atmospheric density 
measurements were not available. Therefore, the modeled descent rate does not match 


precisely with the measured descent rate but the differences appear negligible. 
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Figure 47. Measured vs. Modeled Velocity 
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Initializing the model at the start position of the flight test, the model's ability to 
estimate position of the system was evaluated. Figure 48 shows the model very 


accurately predicts the flight path of the system under the given wind conditions. 
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Figure 48. Measured vs. Modeled Position 


The response of the system to a control input was identified during flight test. For 

a single control input, a glide ratio of 0.4 to 0.5 was found with a time constant of about 
4-5 seconds. The glide ratio was decreased to approximately 0.2 when two simultaneous 
controls were actuated. This response was incorporated into the system model by adding 
the resultant horizontal velocity components. Recall the flight test results demonstrated 
an apparent decrease in descent rate upon control activation resulting from the decreased 
oscillation angles. This reduced vertical velocity is an apparent benefit to the system 
allowing more time in the air to be controlled to the desired impact point. However, the 


implementation of the point mass model does not predict oscillation angle. Therefore, no 
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vertical velocity change was incorporated into the model on control input. The modeled 


response to double and single control inputs is illustrated in Figure 49. 
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Figure 49. Modeled Control Input Response 
As observed in flight test (Figure 50), a glide ratio of approximately 0.2 for two 
simultaneous control inputs and 0.4 to 0.5 for a single control input result. A time 


constant of approximately 3 to 5 seconds is observed. 
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Figure 50. Measured Control Input Response (Two Controls and Single Control) 


The modeled data closely represents the measured system response. 
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XI. SIMULATION 


A. OVERVIEW 


The individual modeling efforts have been presented. Three major components 
are included in the overall system model: 1) Dynamics Model, 2) Sensor Model, and 3) 
Control System Model. In addition, the reference trajectory generator was implemented 
using the same equations of motion utilized in the dynamics model. This trajectory 
generator has the additional benefit of being a tool for calculating the Computed Air 
Release Point. Figure 51 provides a block diagram of the overall control concept. This 
control methodology was implemented in Simulink® (Appendix C) and integrated with 
the sensor models and dynamics model. Single simulation runs can be executed. It is 


desirable to assess to obtain a statistical base for many trials. 
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Figure 51. Control Concept 
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B. INITIALIZATION 


To obtain the statistical base desired, an initialization file was created where 
numerous parameters utilized in the simulation were selected randomly. The 
initialization file then provided the parameters required for numerous simulations. A 
MATLAB? script file called INIT was created to select the desired parameters. When 
executed, the user is prompted to input the number of simulation muns desired and to 
select the "actual winds" to be used. The program then randomly selects the parameters 


identified in Table 4. 


Possible Result 


Planned Wind File Uniform Distribution any valid wind file from 
the same day as the 
selected "actual winds" 


pi enennoroii 

selects which sensor model is used 

distance from desired release point each oe 

PSIDOT Normal Distribution Mean: 1.89 degrees 
ftom Snags — 
RELEASEALT Based on maximum 9500 feet, 20000 feet 


Uniform Distribution 


SASEED Uniform Distribution 0 to 9999 
seeds used in GPS error models 


Table 4. Simulation Parameter 















Another script file, titled CARLO is then executed. This execution file imports 
the parameters in the initialization file, executes the simulation CARP to calculate the 
predicted trajectory and the simulation COPOINTMASS to perform the simulation. The 
results are then stored in individual data files. Sree statistics are calculated by a 


script titled MONTE called from within the CARLO routine and stored in a text file. 
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c RESULTS 


Appendix A contains the initialized states and results for each of the 600 
simulations conducted. These data include the achieved accuracy improvement of the 
controlled system over that of an uncontrolled system. The simulations produced 
excellent results with an accuracy of 210 feet (64 meters), Circular Error Probable (CEP). 
The total average horizontal error was 309 feet (94 meters) with an average of 15 control 
inputs being required. The maximum number of control inputs for all 600 simulations 
was 33. Appendix B presents sample plots from several simulations. 

Figures 52 through 54 present three dimensional plots of several of the simulation 
results. These figures illustrate the initial release points and flight paths to impact with 
the center of the coordinate system (0,0) being the planned impact point. With the 


exception of a few trials, the system guided to the desired target. 
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Figure 52. Simulation Results, 10000 foot Release, 100 trials 
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Figure 54. Simulation Results, 50: 10000 foot Release, 50: 20000 foot Release 
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To assess each of the randomized variables, the statistics were summarized by the 
time from the predicted winds to the actual winds, the GPS errors introduced (state of 
selective availability (ON = included)), and the distance from the predicted CARP to 


simulated release point. Table 5 presents these results. 
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Table 5. Simulation Results 


The predominant factor microg accuracy of the control system was the time 
from the predicted winds used to establish the planned trajectory to that of the time of the 
simulated airdrop. Figure 55 plots the resultant accuracy with the “age” of wind and 
shows that the impact accuracy is greatly reduced when the actual winds are within 


approximately three hours of the predicted winds. 
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Figure 55. Influence of Age of Wind on Accuracy 

Figure 56 presents an example of the magnitude of wind error for the three and 
five hour old winds. For the first set of winds (three hours), the actual winds closely 
match those of the predicted winds (within 15 feet per second). The second example 
(five hours) shows differences in wind speed up to 30-40 feet per second with changes in 
significant changes in direction below 2,000 feet in altitude. The system simply does not 


have enough drive to overcome differences of this magnitude between the planned winds 


and the actual winds. 
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Figure 56. Wind Comparison 
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The effects of GPS navigation errors are not as apparent as that of age of wind. 
Figure 57 presents the accuracy with the status of GPS errors (Selective Availability). 
This plot shows no direct correlation between the state of GPS errors and the resultant 


impact accuracy. 
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Figure 57. Influence of Selective Availability Errors on System Accuracy 


Evaluating the trials that had impact errors less than 100 feet, a slight correlation 
is observed between Selective Availability status and impact errors. For SA Off tnals, 
the resultant impact errors had a smaller variance than that for SA On (Figure 58). This 
shows that the GPS navigation errors do contribute is the overall accuracy of the AGAS 
system but do not have the same degree of impact as the age of wind data. 

The effect of the magnitude of the offset of the actual release position from that of 
the planned release point also does not show a direct correlation with the magnitude of 


the impact accuracy (Figures 58 and 59). 
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Figure 59. Influence of CARP Errors on re Accuracy 
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XII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


This study demonstrated that the AGAS concept is capable of providing a cost 
effective option for precision airdrop. Simulations demonstrate that accuracies of 210 
feet (64 meters), Circular Error Probable (CEP), can be achieved. The development is 
warranted in going into the next stage of development. Additional efforts are needed to 
optimize the control system in an effort to further reduce the amount of stored gas 
required on the operational system. 

The flight test program provided adequate flight dynamic data for the AGAS 
system. The instrumentation system developed for this program proved invaluable. The ~ 
wind measurement techniques employed by use of a "calibration" parachute system 
provided the best possible wind measurement system for this application. Use of the 
measured GPS ground track velocities as the wind estimate was validated to provide 
sufficient wind measurements. 

The MMLE3 parameter estimation algorithms as implemented in MATLAB386 
are not sufficient for predicting the stability derivatives for this parachute. A numerical 
sensitivity to the magnitude of the particular parameters is suspected to be the cause. It is 
unlikely that these techniques, had they provided a solution, would have produced a 
reliable solution due to the non-linearity of the aerodynamics of a parachute system. 

An efficient Monte-Carlo type simulation was developed using a point mass 
model for parachute dynamics, sensor models for GPS and heading information, and a 


Bang-Bang type control system. 
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The point mass concept for system modeling 1s adequate for this evaluation. The 
position and velocity results from the model matched well with the measured flight test 
results under the same wind conditions. However, to optimize the control system, a full 
six degree of freedom model is likely required to capture the proper heading response. A 
complete set of general 6-DOF equations of motion were developed and presented herein. 
However, additional work is needed to either estimate the stability derivatives for the 
parachute system or develop a physics-based aerodynamic model such as that proposed 
by Colin Tory. 

Six hundred simulations were conducted with randomly selected initialization 
parameters. These results demonstrated that the Affordable Guided Airdrop System, as 
described herein, shows strong potential of providing a low-cost alternative for precision 
airdrop. Three critical factors will drive the final design of the AGAS. First, the 
accuracy of the estimated winds when determining the planned traj ectory is the dominant 
factor in the accuracy of the AGAS concept. Winds of up to 6 hours old (as compared to 
the ‘actual’ winds used in the simulation) resulted in large horizontal errors from the 
desired impact points. Secondly, the rotation rate of the parachute system is important. 
Rotation rates with a mean of 1.89 degrees per second and a standard deviation of 1 
degree per second allowed effective control. If the rotation rates of the production 
system are increased from that, sufficient control may not be possible. Finally, the 
number of control inputs required to achieve the desired accuracy is marginal under the 
current control concept. However, no attempts were made to optimize the control 


algorithms for minimum fuel usage. 
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B. RECOMMENDATIONS 


This study demonstrated the feasibility of the AGAS concept to provide a low 
cost alternative for precision airdrop. However, the success of the final design rides on 
the three critical factors presented above. Therefore, the following recommendations are 
offered for follow-on work: 

1. Complete the incorporation of 6-DOF equations of motion. Nonlinear 
parameter estimation techniques would need to be investigated or a physics-based 
approach would be needed. 

2. Fully characterize the performance of the AGAS concept using the G-12 
and/or G-11 parachute systems. The remote controlled activation technique used on the 
C-9 test program should be applied to the G-12 system as soon as G-12 actuators are 
available. 

3. Investigate optimizing the control algorithms for minimum fuel usage. The 
current methodology provides minimum horizontal errors without regard to fuel 
consumption. The fuel consumption for these algorithms is marginal. A technique for 
evaluating multiple control algorithms 1s the application of a ground based guidance 
computer. The navigation data from the system could be downlinked via radio modem to 
the guidance computer, the control algorithms could determine the desired activation of 
the actuator, and these control commands could be uplinked to the test item. Figure 60 


illustrates this concept. 
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Figure 60. Autonomous Guidance Test Concept 
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APPENDIX A. SIMULATION RESULTS 
Error Control 
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OS ian ios 500 err seas pa 
3610 | s000 [3000 | 25 | 9500 | sae | ea7 [sees [CdS 
a7 7 | 0 | 100 | 1000 | 3.0 | 500 | seas [516 | 2370__| 735 | 11 
eft ~Pto00 [100006 [9500 | -e1e7_| 963 | 2150 | aes] 

To fo | o [181] 9500 | 3004 | 667 [778 | 6 | 20 | 
[2000 [2000 [1.7 | $500_|-3677_[ 64 | 2770| 72 | 14 _1 
mt 0 On roa [meson nfpmeseeau [687 | marvin | zanna mn 
A 

7 fo | o | 08 | ss00 | -s677_| es [1026 [169 | 18 | 
7 [0 | -2000__| 2000 [3.3] 9500 | -2677_|_-84_| 2774 |e «dC 
ast 2000] 2000 [7.0 [e600 -3e64 | oa? | sav | ter 
p46 [3 TY «tno | tooo [24 | 9600 | -se77_| 4 | 481] sad) 
a7 [73000-3000 [21 | e500 | cess [516 | 4176 | aes] 
eft | 2000 | 2000 fos | 9500 | 3864 [ear | sea7 | ar] St 
Das [7p ef tooo | 1000 [03 | 9500 | cess | 516 | 2374] see] 
ro fo [| 0 [18 | 9500 | -s864 | 687 | 776 (36 | 18 
spe fo | 1000 [1000 [08 | 9500 | 6187 | ces | 2147 [ea sf 
sap 1000 | 1000-25 | 9600 | -se77_| es [482] ~sSYSC 
a3 fe dt 2000-2000 [24 | 9500 | 6is7 | ves | 05s toast 
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Hours from CARP Y Release Error Error Control 
Pred Wind Offset Offset Alt Predicted CARP | NoControl | Controlled | Inputs 
ss [| 6 | 0 | 3000 | -3000 | 24 | 9500 | e1o7 | oo | ees | 48 | 13 
se [1 | 0 | 1000 | 1000 | 12[ 9500 | -csca | Gey [| 798 | s8 | 13 
BSE, an ee ee ee 
58 | 8 ee ee 
Re a ee Oe i 





Wi NT G@] — 


a 
[7 fo fo { 0 | 19 | 9500 | -6e8s [S16 [3565 | 108s 3 
1s | 9500 | sea | a7 [77a] tez—«d| te it 


~J 


ee a 
ef 7+ 1 «2000-2000 | 2.3 | e500 | -6eas_ [516 | 5a33__| izi0_| 9 
ea_f 7 +1 0 | 3000 | -3000 [0.0 | 9500 | -6eas [sie | 41s | aaa [5 
Te pe | 1 | 1000 { 1000 [23 | 9500 | -e1e7 | 963 | 2ta7_|_300_| 7 
ee [7 | 1 | -2000 [2000 [27 [| 9500 | _-608s_| S16 | 5428 | aaa 10 
er [1 | 0 | -2000 [2000 | 26 | 9500 | -seea_| ear | s2a1_| 26 | 13 
ee [1 | 0 | 3000 | -s000 [2.4 | 9500 | -s8e4 | 687 | 3005 | 35 | 14 _ 
es [7+ 1 [3000 [3000 [1:9 [9500 | -6eas [ie | ates ted Ct 
[7a [7 | 0 | 1000 [ 1000 [25 | 9500 | -ecea [sie | 2374 | 924 | i2__ 
|| tO [ey OL a2 Vies00 i |mm-sae4 yy) ce7 al 7740 |e 
wz psd 000 8000 6 | 9600 | 3677 [4 | 4505 «| Od) 
eo s000 | -s000 5 | 9800 | sees | “ste [ater | sao Jo 
7a f+ 1 ~3000 | 3000 [16 | 9600 | -61e7_| 963 | 4008 | tea__| 12 
cm psf 1 [oO | 0 | 04 | 9600 | -se77 | 84 [ 106 | 03 | 16 
we [6 | 0 | tooo | i000 [16 | 9500 | -e1o7 | 963 | 2ta7_| 14a [ 11 
[7p 3 to [3000 [ -s000 [20 | 9500 | -se77_| 04 | asa | 45 | 15 
A A 
7s [6 ‘| 1 | 3000 2000 [33 | 9500 | -6187_| -963 | 4005 | ser | 14 
eo [7 -+| 0 | 1o00_[ 1000 [15 | 9500 | -ccee [516 | 2378 | s17_| 13 _ 
met] scm. 1 amo" eNO) ae25. ues00 (m-covrmn| 04 luust025 on |emtco. nn oa 
pez pf 6 | 0 | 1000 | 1000 [26 | 9500 | -e1s7_| ces | 21sa | se | 9 
asf 7 *dY~—t~C:«C 2000 ~[ 2000 [28 | 9500 | -coae [516 | sama | 933+ = 
[To [0 | 09 | 600 | -seea [eer [775 | 90 ~+| 16 
[sooo [3000 [2.1 | 9600 | cess [516 | 4176 | ee | om 
a 
et free Oni = Snap PO.8 Tie G00 “1p > S677) 84a] megO2 mone 20S ater 
css 1 tf 2000 [0004.1 | 9500 | 3864 | a7 [324s —its2_sd ts 
[SS oe | eee San [oe Oommen] ame | 920 6 tu 500 we [w= C077. | 204-0] cer 1024 Oe | eee cee | ta 
A 
Sons ater ees ar tes 
a [1 {0 | 1000 | 1000 | 79 | s800 | -s86a [ear [80s | —a7_—=id;s=Ci‘O 
[3 | 0 | 3000 | -s000 | 1.2 | 600 [| -ser7 | a4 [4600 [| © | 16 
4 [7-1 | 3000 | 3000 | 7s | 500 | ees [S16 | 4i7i [ees 
[6 | 0 | 3000 | -s000 [28 | s600 [e1e7 | ses | aese | ies) 11 
5 | eS Epa 1 oe [ ew One NOM MOTH [RSS | NES0774[™-G4 [10247 [1688 | aa 
7 [6 | 0 | 3000 | 3000 | 26 | S600 | -ete7_| ses [__4ee7_ [813 
6 || aT | COC S00 [= 7 7a [=A 1024 | e222 
<a 
oO | 1 | 1 | 3000 | 3000 | 7.6 | 9500 | -seea_| 687 | 3003 | 333 | 11 
107 
103 a o 
I 


~J 


~] 


OO} © 
mil 


i 1_|_ 1000 [000 9500 | Tae |-867_| 588 [ase | 7 
105 1 0 1000 1000 | 9500 = 4. -=7283 -867 214 18 












0st 000 000-800} ares at [aise ts 8 
A SO = SZ 
Pos [2 | 1 {| 0 | 0 | 08 | sso | res [or | oo | 23 | 13 | 
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Error Error Control 


Hours from CARP X | CARP Y Release 
Pred Wind ib ae Alt Predicted CARP | NoContro! | Controlled | Inputs 


i a a Tos | 9500 | rae | a1 | 573 | 48 | 4 
a aC ao 
Ciat_[ 1 | 1 [3000 | -3000 [07 |s600_| -7283_| 807 [_a1az__[ _252—=«d|s=Cit 
riz] 5 | 1 | 1000 | 1000 | 08 | 9500 | aaa [ a1 | atte | a7 | 8 
Tis | 1 | 1] _3000_{_-3000_| 3.3 [9500 | -7283_[ 867 [| aa1_—[—sse—idYSCi 
a LT A DZ 
Soon oe a ma Oe) Ome 06) S500) -a74sn|. sie 2008 | ces ee 
Tie} 1 | 1 | 3000_|_-3000_[ 1.8 | 9500_| 7283 | -ee7 | 4145 =| =r] SO | 
[air [5 | 1 | 3000 { -s000 [07 | 9500 | -47as [a1 | 6723 | ees | 3 
[a yA = SS OC 
TN LL B= (S 
Oo So 2 
Paz i om, oF | 0 | 20 \menonesomy elle eee. |. ie nn 
Bie. 2 | _o) |p -awo | soon 2s [eo | reas | ot) ses, |e 
ee oe a, ca) 2000) 000 | 070 | 500 | -araon | at BE ICE | oS 
ene oe Oe fe Oe ft Oeerle Lar [peS000n | -recam) OT ai ees eee te 
| 2 | 0 | 1000 | 1000 [76] 9500 | -75a3 | -o1 | io | 169 | 15 +4 
[2 | 1 | 1000 | 1000 | 06 | 9500 | -7sss | 01 | 1365] 98 | 10 | 
7 | 5 | 0 | -2000 | 2000 [05 | 500 | 47a | 41 | 2103 | 27 | 11 
en) wien | om) On) ore |) ar jrreccon lp -raes| 887 ln OD le 
[wt 5 | 0 | 0 | 0 |19| 900 | 4743 [| a | 1908 [ «1 [| 10 
liso | 5 | 1 | -2000 | 2000 [12 | 9500 | 474s] 41 | 2192 | 246 | 12 
131 20 
cea |e el) | Ome OW) 2A: |mCcOON IN -7593N|) Simla S| em 
css | sof see | -anoo [os | ss00 | ares | at | rez rors 
ise [eo sao | -a0oo | 32 | ss00 | ress | at oeat [te [8 










ee “20mPG)) 0 z1 | 9500 | -7283 | = Ee ee 
136 588 | 2000 | "9500 | -4748 | 41 2188 237 18 | 









raat [Sf a00 at] ss00 | ares [ a tas Pes 8 
Tise_[_2 | 0 _| 3000 | -2000-| 1.2 [9500 [7583 |_-o1_[3as0_ [134 [10 
rise [5 | 0 | 10007 | 1000 [0.8 | 9800_[-araa_[ai_[ are | 2—«d 0 
rao [2 | o [-o |_o | 18 | 900 | -ree3 [or [ota 
ria [ 2 fo | _o | 0] +5 [9500 [7803 [91 [ 90247902 
Diaz [ 2 | 0 | 2000 [2000 [1.0 _[s600_|-7533_ [81 367a [tar 
Tas_[___1 | _1_|_-2000 | 2000 [2.3 | 9500 [7283 | ser [3186 [a2 
[14a_[ 5 | 1 | -2000~| 2000 [33 [9500 [-araa_[a1_[aie3 [as te 
Caaf 5 Et 20002000 [43 | 9500 -aras_ [a1 2te2 [0a 8 
riae[ tf _t000 [0004.0 [9500 _|-7283_[ 867 [02 [20 
pas_[ 2 [0 | 3000_[-3000_ [2.5 [9500 | -7503[-e1_ [sass [tae 
Dis[ 2 | 0 | 3000 | -2000_[1.0-[~9500_|-7e93 [a1 [saan tt iT 8 
Piso fs | 4 | 0 | 0 | 06 | 9500 | --ara_[~ai_[ tena [ sts ST 
Disa f 5 0 | so00 | -s000 [+8 | 9500 | are [a1 | s72a [es 
css [1 | 0 | 1000 | 1000 | 23 | 9500_[ 728s | 867 [574 [tats 
risa [ 5 [1 | 3000 | -s000 [4.3 | 9500 | -a7aa_ [a1 [5725 [toes 
ras5_[ tt «1000 [0004-7 |9500_[-7283_ [867 [S80 [es a 
p1s6| sof tooo | F000 |. | 800 ares at [azo aoe 
is7_[ 1 [1 20002000 4.8 | 9500 _| 7283 [se7_[ sie fs 
Piso [5 | 0 | 10001000 [4.2 | 9600 [Aras [ai ano [es 
Creo [5 [1 | 0 | 0 | 05 | 9500 | -ara_[—a1_[ 1903 [32s 
Pier | 5 | 1 | sooo 3000-23 | 9500 [aaa [a1 [7a at 
[iez [2 | 0 | 1000 [1000 | 0.8 | _9500_[-7588_[ or [tava 
ies [1+ 1 | 2000 [2000 [49 | 9500 [7283 [867 [aie iar 
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Error Error Control 


Hours from CARP X | CARP Y as 
Pred Wind Offset Offset Predicted CARP | No Control | Controlled | Inputs 


a 
Liha ae [24 | 9600 | -72e3 [867 | se1_—«[| ~—Sto2—Ssd|SCt CY 
a 
| 167 | 2 | 1 OE 3000 =| -300G 1.6 e500 1h e75saRR tl) * -91 i 3432) SS ee 
ree [i] tC «2000 =| 2000 | -02 | 9500 [| -7283 | -ec7 [ 3132 [| a0 CTC CC 
pies [7 ~[ 0 [ 1000 | 1000 | 20 | 9800 | -72e3 | -ae7 [507 [| 279 ~=«|,sist]2 
Piot 5 | 1 | -2000 | 2000 | 34 | 9500 | 4748 | 41 | 2168 | 29 [| 18 
Pitt 2 [| 1 | -2000 [2000 [29 [ 500 | -76s3 | -51 | ses [77 =~ ~SOt4 
Pi72[ 2 | 1 | 100 | 1000 | 03 [ 9500 | -7603 | 61 | 1355 | 104 | 11 | 
pis {1 | 61 ~<(| 2000 | 2000 | 18 | e600 [ -72as | -ee7 [ 3ic6 | tea TCS 
Piva 1 [| 0 |{ 1000 | 1000 [12 | 9500 [ -7263 | -a67 | 582 [| ter | 7 
is [5 [ 0 | 1000 | 1000 [os | 9500 | 4748 [41 | 3203 —| So aoSTtC* 
Pie { 5 | 1 | 3000 | -3000 | 06 | 9500 | 474s | 41 [ 5724 [ 1129 [ 4 
pia [2 {| 0 [ 1000 | 1000 [21 | 600 | -7es3_ [or [tsze ta Tat 
Pie [2 [8 [ -2000 | 2000 14 | 9600 | -7es3_ | 1 [3672 [ter 
Pit 2 | o | 3000 | -3000 [719 | S600 | -7533 | -81 | 3433 [148 | 17 
Pteo [1 | 1 | 1000 | 1000 [15 | 9800 | -72es [sey [500 [38312 
ter f[ iT 1 =< 3000 | 3000 | 14 [ 9500 | -72e3 | 867 | ai41 [sea [12 
rez tt ~~] «1 «=| «=3000_:~+|| 3000 | 22 | 9600 | -7283 | ser [| 4747 | 194 | 16 
is [5 [ 0 | 3000 | -3000 [26 [ 9500 | ava [a1 [ S722 | 7er_ | 9 | 
peat SCT «St =| ~-2000 + 2000 | 40 | 9500 | 4748 | 41 [ 2io1 [| 276 [ 19 
wes [SCOT St ~=*[«~-2000~*[ 2000 | 05 | 9600 | -474a [ 41 [ 217 [ 315 | 12 | 
pias] 1 ~—| «61 ~<[~-2000 *| 2000 | os | 9500 | -7263 | eer | 3106 | 328 | 10 
Pta7 tt ~CdT St ~=[ «2000 :*| 2000 | 29 | 9600 | -7283 | eer [3106 ~| 183 SCs| 
psf td] OT (CO CT (OC 45 | 9500 | -7283 | 867 | 872 | 265 ~+| 13 
SSR So en eee Ee eo eee 
Piso [1 ~{ ~0 ~| ~-2000 ~| 2000 | 22 | 9500 | -72e3 [ ser | siot | 24 | 73 | 
Pit [tT tf 2000 | 2000 [3.4 [9500 | -72e3 [ -8e7 [31390786 
[iz] 2. | a7 | 0 | 0 | 26 | 500 | -r5o | or) 608 | 17? (on tGiul 
piss f 2 | o {| 3000 | -3000 [38 | o500 | -7633 | or [3435 [| i7 | 17 
Ciee | 2 dT 1 =~ -2000 ~-[ 2000 | 24 | 9600 [ -7533 [| -91 [| sere | 164 [| 13 
Li a Se a ee ne eee PS ae ee 
(a cr a | | Pee eg Ee RR 
[oia| ae tee fe, OMe” O- i _O” ilraO.7 queeGoee fe -72es sige? smear ar leon ail aay 
(Rd fee a | 
(7 ee ee ee ee 
200 a a Sk ee EL ee 
or [so 2000 | 000 [14 | 20000 | 73236 [1503 | io8es iva [aT 
202 0 il On |r GOmNNIEO.2 |lze0eg: [maszaamintees [mee S0625 | 5 25071 Joi 
et 7 | 1 ~[ «(3000 [3000 | 20 | 20000 | 14857 | 1019 [ass [274 
2045 [Sum] 1 = |e 3000 [i 3000 | 2.115 | 20000) [16319-1584 [7 769 | e363 fe 2a 
25 fs fo | 3000 [3000 [a7 | 20000 | i7o12 | -2036 [zen 718 
P2003 Cd «OSC «2000 :~| 2000 [30 | 20000 | i7012 [2036 [ 5e20 [10 |e 
|_207 [1 [000 | 1000 [3.0 | 20000 “| “rvar2_ | -20ae | 31e5 | 32 30 
206 [St 20002000 [| a0000 | tests | 1584 | esa7 ait te 
20 [5 | 0 { -2000 | 2000 [22 | 20000 [13238 | 1503 [~ foeea [156 [7s 
[20 | 7 | 1 | 1000 | 1000 | 20 | 20000 | 14857 | 1019 | e7ei_ | 544 | 20 
fen! 3 1 [pee | 1000 sligt-za|200e0, | tzo12)=20c6 fa SiO ta |S ae 
| 212 Fo | 3000 | -3000 | 26 | 20000 | 14857 | 1019 | 2388 | 205 | 23 
| 





TS 


On 


213 [ 5 | 1 | 3000 | -3000 | 13 | 20000 | 16319 | -15e4| 737 | 406 | 20 | 


214 | & ee 1000 1000 a 20000 | EIASS a P 101Ie | C783 5 | 7 33 
| 215 7 3000 -3000 20000 | 14857 ; 1019 | 2336 226 21 

















~] 


Ss 
a7 {9 [1 | -2000- | 2000 [29 | 20000 | 13238 | 1503 | 108es_| 1320 | 2 
Cais | 9 | 0 | 1000 | tooo [1.9 | 20000 | 13238 [1603 [8113 | 602_| 17 
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Hours from Error Error Control 


CARP X | CARP Y Release 
Pred Wind Offset Offset Alt Predicted CARP ; NoControl | Controlled Inputs 


ras] 3 0 | 3000 | -s000 | 18 | 20000 | 17012 | -2036] 1268 | 35 | 20 | 
ra [7 | 0 | 1000 | 1000 [4.0 | 20000 | 14857 | 1019 | e777 | 130 | 16 
Tai [ 7 | 0 | -2000 [2000 [21 [ 20000 | 14857 [1019 | eses_—| 73 =~dY~Ctes 
pam [3 | 0 | 3000 | -s000 | 26 | 20000 | 17012 | -20s6| 1261. | 121+) ~30 ~~ 
[aa 7 | 1 | 0 | 0 | 35 20000 | 14867 [1019 | escs | 18 | 33 ~~ 
[a [5 | 0 | 3000 | 3000 | 12 | 20000 | 16319 |-1504| 738 | se2_~| 7 
[as [3 | 1 | 3000 | -3000 | 27 | 20000 | 13238 | 1603 | 3651 | 76 | 24 
13 
Tze [—7_| 0 | 3000 | -3000-[~0.4 | 20000 | 14857 | 1019 | 2333 | 189 | 10 
19 
LS aa eS 
= ane] Ome] 0 ozs "ma000 WI teston] Sedum 3726p | ae al ee 
Peete — | One|, Oe] 0] 27 alezan00a| igsemelaota |) ses |e | aes 
[za [3 | 0 | -2000 | ~2000 [26 | 20000 | 17012 | 2036 | sel7_ |. 106 | 2 
235 sf ro00 0003.1 | “20000 | “arora 2086 [sie [70 27 
[236 | 7 {| 1 | -2000_| 2000 [25 | 20000 | 14857 | to19 | ses izes 
27 [3 | 0 | -2000 [2000 [3.0 [20000 | 17or2 | 2036 | _sez1_| tsi_+i| =o 
me | 5 | 1 [0 |~0 | 06 | 20000 | 16319 | -1984| 3722 | 477 | 19 
eae | ae| Owe] Ome) oiSa(a20000 9] 170lam| 2036) Abe; | am | mw 
P20 [ 3 «df 1 | oo | t7 | 20000 | i7or2 | 2086] 2800 | ~—ossdYSCt 
Peat fs | 0 | 1000 [1000 [0.2 | 20000 | teste | -15e4 | sa03[ 22 | 19 
Paap 7 fo | o | 0 | 21 | 20000 | 14857 | 1019 | 0504_[ 92 | 28 1 
Tas fs ~ o | 0 | 0 | 08 | 20000 | 13236 | 1503 | 6055 | 7oe | 68 
24a [7 | 1 | -2000_[~2000 [3.3 [20000 [ 14857 | tore [9388 | s02__| 28 
es |e [tooo [ooo [1.5 | zo0c0 | rvore_| 20a | sz00 [235 [20 
! 2.4 20000 14857 1019 6566 
aarp oo oF oo sie ise | 3808 | a0} 
[24s [5 | 0 | 1000 | 1000 | 09 | 20000 | 16319 | -1584| 3606 | 480 | 20 
[2s [5 | 0 | -2000 | 2000 | -0.3 | 20000 | 16319 | 1584 essa | _242*| 10 
[20 [7 | 0 | 0 | 0 | 148 | 20000 | 14857 | 101s | ses | 135 | 19 
[21 fs | 0 | 1000 _| i000 [0.8 | 20000 | 16319 | -1584 | 806 | 355 2 
Dae [9 | 1 | 1000 | 1000 [25 | 20000 | 13238 | 1603 | 6113_| 673 | 21 
[26 [9 | 0 | 1000 | 1000 | 12 | 20000 | 13236 | 1603 | e111 | 803_| 12 
[28 [7 | 1 | 3000 | -3000 [08 | 20000 | 14867 | io1a | 2336] te7__| 12 
[20 [9 | 1 | 3000 | -3000_| 1.8 | 20000 | 13238 | 1603 | s626| 203 ~«| +18 
Paez [7 | 1 | 3000 | -3000 [15 | 20000 | 14867 | 01s | 23ar_| 53 | 
2 a 
T2eaf9 | 1 | 1000 [1000] ~o0 | 20000 | 13236 | 1503 | 8i13_| aa] 7 
as [00] a0} S000 | estore | ass 
[2s [5 | 1 | 1000 {1000 | 27 | 20000 | 16319 | -1564| 3807 | 261 | 25 
[27 [5 | 0 | 1000 [1000 | 26 | 20000 | 16319 | 1584 | 3805 | a72_—+| 8 
Tae5 [3+ 1 ~| 3000_| -3000 | 24 | 20000 | 17012 | 2036 | 127 | 10 | 2 
[269 [7 ~~+| 0 | 1000 | 1000 [18 | 20000 | 14867 | 1019 | ear] 45 | 
rao [3 +t 0 | 1000 | 1000 | 25 | 20000 | 17012 | -20se[stss_ | 51 | 2 
Dari fe | 1 | 1000__[ 1000 |~2.4 | 20000_| 13238 [1503] ei13_—~| sees 
Taz [3 tt 3000-3000 | 1.6 | 20000 | 17or2 | -20se| 1284 [am if 0 
Pars fs«d at -2000-[ 2000 [15 | 20000 | 13236 | 1503 | to8es__ | ea 2 
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Hours from CARP X | CARP Y pat | rasa Error Error Control 
Pred Wind Offset Offset Alt Predicted CARP | NoControl | Controlled | Inputs 
(37a [5 «| 0 ~~ 3000 ~| -s000 | 16 | 20000 | 1es19 | -1684| 738 | 274 =| 22 
[27s [9 | 0 | 3000 | -3000 | 04 { 20000 | 13238 | 1503 | 3623 | 172 | 11 
part 3 | o | 1000 | 1000 | 20 | 20000 | 17012 | -2036| sim | 169 =| 20 
Pave [3 | 0 | -2000 | 2000 | 22 | 20000 | 17012 | -2006| sere | 110 | 18 | 
Pare [3 | 1 | -2000 | 2000 | 33 | 20000 | 17012 | -20se| ses | 71 | 25 
P20 [9 | o | 1000 | 1000 | 15 | 20000 | 13238 | 1603 | e112 | 75 | 14 ~*| 
Zar | 5 {| 0 | -2000 | 2000 | 25 | 20000 | 1es19 | -1584| 6540 | 176 ~| 21 
P22 {3 {| 0 | 1000 | 1000 | 26 | 20000 | 17012 |-a0se| sia | 58 | 30 
Pass | 7 ~—(«d|~«( «1000 ~| +1000 | 28 | 20000 | 14857 | 1019 | e7e2 | 110 .| 2 
Pasay 9 ~~«| «~~ ~3000 +| -s000 | O98 | 20000 | 10208 | 1603 | 3825 | 269 | 11 | 
25 | SCT Of «1000 ~| 1000 | 21 | 20000 | 10208 | io03 | 8113 | 775 | 18 
ase | 7 -| 0 | 1000 | 1000 [21 | 20000 | 14857 | 1019 | 67ei_| 171 | 26 | 
[27 [S| «Of ~3000 ~| -s000 | 08 | 20000 | tes19 | -isea| 738 | so | 14 | 
288 [7 ~+| 0 | -2000 | 2000 | 11 | 20000 | 14857 | 1019 |  ooae | 1652 ~| 14 
Paes [3 {1 | 2000 | 2000 [29 | 20000 | 17012 | -20s6| 5820 | 32 ~+| 23 
Pao ft 3 {| o | o | o | 17{ 20000 | 17012 | -2036/ 300 | 107 | 2 
Parts «dT ~~ ~—3000_~| -3000 | 22 | 20000 | 17012 | 26, i283 | 69 | 2 
P22 ft 9 | o | 3000 | -2000 | 19 [ 20000 | 10238 | 1503 | 3827 | 26 | 19 ‘| 
rast 3s YT Oo CT CO SCT =~ SCdTst2 | 20000 | 17012 | -2006| 2006 [| 2 | 20 | 
P24 ft 5 | 0 | 3000 | -a000 | 13 | 20000 | tesa | -1564[ 735 | 182 | 20 
P25 [3 | o | -2000 | 2000 | 23 | 20000 | 17012 | -20s6| sez | 74 | 2 
Ta [3 | 1 | 1000 | 1000 | 09 | 20000 | 17012 | 2096, 31% | 253. | 15 
P27 [7 | 1 | 1000 | 1000 | 26 | 20000 | 14857 | 1019 | 67a | 34 | 2 | 
Pa [3 | 0 | 3000 | -s000 | 24 | 20000 | 17012 | -2036| 1288 | 72 | 30 | 
Pa [9 | 0 | 1000 | 1000 | 1.0 | 20000 | 13238 | 1503 | 8110 | 637 ~| +10 
P30} 9 | 1 | -2000 | 2000 [16 | 20000 | ta238 | 1603 | toees | 1117 | 12 
zor | 7 {| 1 | -2000 | 2000 | 25 | 9500 | cose | 516 | 548 | 1200 | 9 
P30; 7 | © | -2000 | 2000 | 12 | 9500 | cess | sie | Sar | 1193 | 4 
Ps [7 | © | -2000 | 2000 | 24 | 9600 | -seea | ea7 | s28 | 4 | 15 | 
F304 [1 {| 0 | 1000 | 1000 | o8 | 9600 | -sees | cer | veo | 31 | 14 ~* 
P30 [6 | 1 | 1000 | 1000 | 19 | 9600 | 6197 | 063 | 2150 | «96 | 10 
ir a ee ee ees eee 





GW] WO] OI NT On 


WIN] GO} G)] Oi] © 


[7+ |_1_[ 1000 [T0014 | 3500_| 3864 [ea7 | 600] e2 ~| 14 
S08 00 0003500] 367 | as a 
S005 |e Se Ome 0 w|n5_ |mes00 | sora] Sedma 10255 | 628 | noe 


G) 


Pacem 
VEE ee ee ee ee ee 
a7 [1 | 0 | 3000 | -s000 [ 12 | 9600 | -seoa | ce7 | coer | 4 | 6 | 
3127 | a | ONO | ONO: NIMCCOO!e|™-come| eam 10255 (OR) a al 
A Sa RY A, MI le 
(34 [6 | 7 | 3000 | -s000 | 26 | 9500 | 6197 | 003 | 4605 | 310 | 12 ~~*| 
315 no} 3500 ase] 316} sas | 05 : 
se | 6 | 1 | 3000 { -so00 | 22 | 9600 | -61o7 | ses | 46605 [| 137 _~| «14 
rai | 6 | 1 | 1000 | 1000 | 28 | 9500 | e197 | oo | 2151 | 262 ~| 11 
| 316 10 
ss [6 | 1 | 3000 | -s000 | 12 | 9600 | 6107 | 003 | 4667 [| 690 | 6 
ini 
srt 3 | o | 3000 | -s000 | 08 | 9500 | -sor7 | 84 [| 406 | 90 | 10 | 
[F222] SO || OR 1000 eli TapOmn mate NcoOOry|mecacamli coz en yocae |e uz |e 4 
(3s |i ~SC| «1 ~C«L:«C2000:*([ 2000 + 34 [9500 | seca | cer | s2o | 205 =| 
324 4 3000 Sa 9500 | -3864 | 687 ee 11 


1 
325 3 PO ll = 3000— a4 ~ 9500 | -3677 -84 4590 16 
336-000 220 | a5 an a 
sr [1 | 77] 0 | o | 149 | 9500 | -sees | 687 | ve | 33 | 16 | 


OO Da =z eS 
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Error Error Control 


Hours from CARP X | CARP Y ne 
Pred Wind a aan Predicted CARP | No Control | Controlled | Inputs 


a 24 | e500 | 608s [S16 | 3664 | ree |S 
(a es ce ee ee 
rss [6 | 0 | 2000 | 2000 | 23 | 9500 | 6187 | oes | 5053 | 1050 | 8 | 
Dss2 [3 | 1 | 1000_| 1000 | 25 | 0500 | -se77 | -e4 | 485 | 198 | 18 | 
a TN (= OO 
saa [ 7 | 0 | 3000 | 3000 [2.7 | 9500 | -eeas [ste | _ai7a—-| ava —id)sCO + 
[33 [6 {0 | 3000 | -s000 | 42 | 9500 | ere7 | 003 | 4670 | 313+] 17 
rss [3 [0 | 0 | 0 | 24| 0500 | se77 | -e4 | 102 | 67 | 18 
P37 [6 | 1 | 1000 | 1000 [18 | 9500 | ster | 63 | 2162 [640 -~| 10 
3s [7 | 0 | 3000 | 3000 | 28 | 9500 | -seae | ste | “181 | 505 | 11 
3s [6 | 0 | 1000 | 1000 [1.9 | e600 | -ere7 | 003 | 2te2 | vee | 9 
psa_[ a oor | aso fase ea? |r at 6 
Tsar [ 6 | 1] 0 | 0 | 23 | 900 | oer | os sa | oa | 
a 
i 
[aaa [ttt | 3000 | 3000] 3.1_-|9500_| see | ca7 [| see7[te0~d|~CitS 
[345 [7 | 0 | -2000_ [2000 [37 | e500 | oes [S16 | saa0_| 1156 | 12 
rae [7 fo [oOo | 0 {13} 200 | cose [ore | 3563 | 870 | 1 
[34s [ 7 | 0 | 2000 | 2000 [1.0 | 9500 | 088 | 816 | saze_| izes 4 
Tae [ 7 [1 {~~o | 0 | 25 | 9500 | oss [S16 | 3662 854. | 10 
350 [7 | 1 | 2000 | 2000 | -05 | e600 | -ecas | 816 | saz ties | 3 
rst [ 7 | 1] 0 | 0 | 28 | e500 | cos [S16 | 3563 | 648 | 12 
Pass [1 | 0 | 2000 | 2000 [2.1 | 9500 | -se64 | c87 | sae | 42+ 13 
[35 [1 | 0 | 2000 | 2000 [16 | 9500 | -se6a | e87 | soe | 32. -i| 14 
Psst 1 | 1] 0 | 0 |17| 900 | sees [car [79 | 124 | 7 
Dss7 [6 | 1 | 2000 | 2000 | 06 | e600 | -ets7 | 003 | 5055 | 14a [1 
[ae [ 1 | 0 | 3000 | -s000 | 0.2 | 9600_| -seea | ca7 | s0e0 | 69 | 12 
Tso [3 | 0 | o | 0 -| 21 | 900 | -so77 | 64 | 1028 | 52 | 16 
Psst] 3] 0 | 1000 | 1000 [13 | 9600 | -se77 | ea | 40 [34 | 7 
sez [7 | 0 | 2000 | 2000} 05 | 9600| seas | S16 | 542 | te00_ | 1 
[ses] 1 | 0 | 1000 | 1000 {29 | 9600 | -seea | ca7 | 793 | 50 | 14 
Ps [6 | 1 | 3000 | -s000 | 18 | e600 | -c1e7 | -ocs | acco 572 | 1 
Ps [1 | 0, o {| 0 |18| 9500 | see | ca7 | 77 | 36 | 18 
[ser [3 -+| 0 | 2000 | 2000 | 16 | 9500 | -s677 | 64 | aes | 57 | 12 
SE | 
Ps [1 ‘| 0 | 2000 | 2000 | 1.6 | 9500 | -seea | ea7 [| soa | 46 | 14 
Dara [ 6 {| 1 | 3000 | -s000 | 02 [S600 | rer | 963 | 471] 7s] 4 
Pam [ 7 | 0 [0 | 0 | 08 | 900 | esas [S16 [3565 | 68e_| 12 
raz [6 | 0 | 1000 | 1000 [21 | 9600 | 61e7 | 063 [2148 | 400 ~+| 9 
Das [3 | 1 | 1000 | 1000 | 06 | 9500 | -s677 | -e4 | 462 | 30 | 2 
[ara [1 | 0 | 1000 |~t000_| 1:5 | 9500 | -se6a | 0a7 [vos 19 ~+| 15 _ 
[a7 [7 | 0 | 3000 | 3000 | os | e500 | -ecas | S16] aie2_ | 440 ~+| 12 
Dare] 3 | 0 | 1000 | 1000 | 07 | 0500 | -se77 | es | 481 | 72 | 16 
[a7 [6 | 1 | 0 |~0 | 08 | 900 | eter | oes [3470 [a7 | 6 
Vara 7 | 0 | 0 | 0 | 14 | 9600 [esas | ste | seer | sa3_ | 8 
rsa 7 [0 | 1000 | 1000 [2.0 | e500 | -eeas | S16 | 2372 | 700 +| 8 
Psst [6] 0 | 1000 | 1000 | 29 | e600 | e1e7 | oes] 2163 | 271 | 10 
Tse [1 | 0 | 2000 | 2000 | 1.6 | 0600 | -saea | ca7_| 3243] 40 ~*+| 13 
Psst 6 | 0 | 3000 | -s000 [03 | 9500 | -e1e7 | cos | asco | 206 | 6 





=] WW! Gd) — 


me 











Hours from CARP X | CARP Y Release Error Control 
Pred Wind — ae Alt Predicted CARP | NoControl | Controlled | Inputs 
La (19 | 9500 | -see4 | a7 | 775 | o7 | 18 
ee 
see [3 | 1 | ‘000 | 1000 | 22 | 9500 | -se77 | 64 | 484. | 121 ~| 16 | 
37 [1 | 1 | 3000 [ -s000 | 23 | 9500 | -seea | 87 | 3061 | 194 | 14 | 
386 [3 | 0 | -2000 | 2000 | 31 | 9500 | -se77 | 64 | 2775 | 27 | 7 
ses [3 C«dT ~*~ 2000 +| 2000 | 13 | 9500 | -se77 | 84 | 2772 | or | 16 | 
30 [1 | 1 | 1000 | 1000 | 21 | 9500 | -s8ed | cay | 7or | o2 | 13 | 
rai [7 {1 | 0 | 0 | 18 | 900 | cose | 816 | 3564 | 10 | 14 | 
Pszt 3s | 1 i| o | 0 | 13 | 900 | -367 | 84 | 10 | 321 | 14 
Paes ft Cd] 7 [| OO «| «40 | +9500 | sees | cay | 779 | 198 ~| 10 
rsea [3 ~i| o | -2000 | 2000 | 16 | 9500 | -se77 | -@4 | 27e9 | 24 | 11 
[3a [1 | 1 | 3000 | -3000 | 04 | 9500 | -see4 | c87 | s0es | 12 | 8 
rsa [7 | 0 | -2000 | 2000 | 22 | 9500 | -eess | sie | 5429 | 975 | 7 
see [3 | o | -2000 | 2000 | 07 | 9500 | -se77 | 84 | er | 4 | 14 | 
Paes} 6 | o | o | 0 | 25| 900 | e197 | 063 | 3472 | 542 | 11 | 
aoa ft 7 | 1 | 0 | 0 | 31 | 9500 | esas | 516 | 564 | 85 | 16 
Tao [5 | 1 | 3000 | -so00 | 10 | 9500 | -474a | a1 | 57a | 933 | 4 
(a2 [5 | 0 | -2000 | 2000 | 20 | e500 | 4748 | 41 | 210 | 233 | 13 
Pas [1 | 0 | oO | 0 |17/| 000 | 723 | se7| 873 | 212 ~| 10 
[aoa [5 | 1 (‘| 1000 | 1000 | 07 | 9500 | 4748 | “1 | 3228 | 324 | 5 
[aos [1 | 1 | 3000_| -s000 [30 | 9500 | -72e3 | -e67 | 4138 | 364 | 14 
Paoe [S| St ~*~ 3000_~| 3000 | 78 | 9500 | 4748 | 41 | 5723 | 700 | 6 
Tar [5 | 0 | 1000 | tooo | 24 | 9500 | -a7as | 41 | 3208 | 369 | 15 | 
Taos {7 | 7 | 3000 | -s000 | 17 | 9500 | -/283 | -ee7 | 4198 | 322 | 10 
Paos {1 | 7% | 1000 | 1000 | 12 | 9500 | -72e3 | ser | 582 | 62 | W 
Pato 2 | 0 | 3000 | -s000 | 22 | 900 | -reas | -o1 | 343 | te7 | 14 
Car] 2 [| 7+ | 1000 | 1000 | 20 | 900 | -7533 | -o1 | 1346 | 123 | 14 ~~ 
Taiz {2 | 0 | 1000 | 1000 | 37 | 9500 | -75a3 | 01 | 1956 | 153 + 13 
ROTO 0 00 rs er a a 
Tara [2 | 0 | 1000 | 1000 | 41 | 9500 | -75a3 | -o1 | i331 | 143 ~~ 20 
Cais [7 | 7 | 3000 | -s000 | 34 | 9500 | -r2es | -ae7 | 4137 | 374 | 13 | 
pais, | SS |e Lgm[ees000™=|p 2000 ut |e os00"s| weamcom| le comes eS | 
a7 | 5 [0 np tO00 fs 1000 O2 essOO hp -a748 41 aaio | 367) 
4195 [20 | 1 eal] 2000 HW) 2000TNINNT.S 9] "9500. )m-7539 1m] 0-01 NJmriee71_ || 221 | Onn 
ras [5 | 1 | 3000 | -3000 | 17 | 9500 | 4748 | 41 | 5724 | 1004 | 6 
Tao ft 5 | 1 | 1000 | 1000 | 19 | 9500 | ave { 41 | sate | 409 | 6 
21 i Bd Poe ed 2 
fg 42200| as 2 a | | DOO pe =9000 | | NB @SO0 | NE 75S40N [N-OT a) N=N29 me | D250 | mae 
Paes | 2 | 1 | -2000 | 2000 | 11 | 9600 | -rea3 | -o1 | se77 | 208 | 13 | 
[aoa tT 2. | 1 | -2000 | 2000 | 16 | 9500 | -75a3 | -o1 | 3671 | 109 | 12 ~~ 
ras [1 | 1 | 3000 | -3000 | 1@ | 9500 | -rzes | ey | ai | 82 | 11 | 
Pas [5 | 1 | 1000 | 1000 | 21 | 9500 | 4748 | 41 | s196 | 195 | 18 
[aarti | 1 | 1000 | 1000 | 29 | 9500 | -72e3 | -ee7 | 566 | 124 | 20 
(Za ed ed ee ee 
Taaef[2 [| 9 | -2000 | 2000 | 01 | 9600 | -rea3 | -o1 | sere | 140 | 13 | 
paso [2 ~+| 0 | 1000 | 1000 | 25 | 9500 | -7ea3 | -o1 | 1305 | 166 | 13 | 
(em 2 lp 1 1000 weiGnomiay uaccd0 =jm recom nomi aac | S40 NG 
4325 | 20 | et | OOD 1000 | 2.9 u|NROS00.n [meTess B[eeOTE| w 1954 [5 2510 | os 
ass [oo] 82 | een] rzes _| seer [eee [ss 


-2000 | 2000 = 9500 | -4748 | 41 | 2186 | 130 15 
-2000 | 2000 9500 | -4748 | 41 12 





8 
St 




















rie | 319° 


Loco {ao 7 sc ras ser | 8 
[0 _[- 1000 [1000 | 0.6 |__s500_| 7533-81 | 1352 | 141] ~15 | 
ass [5+ 1 2000 [2000 [23 | 9500 | —47ae_ [a1 [| 2190 | 16 | 13 





| = 
436 



















EEEE 
W 
~] 
tk 


un 
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Hours from enonl SAREE CARP Y ape Error Error Control 
Pred Wind Offset Offset Predicted CARP | No Control | Controlled | Inputs 
a 
| a P13 [ 9500 [ -7533 | -o1 | ee | 128 | 13 
<5 ec a a ea 
Paaz tt | 1 ~~ -2000 ~[ 2000 | 08 | 9500 | -72e3 | Gor | aise | it7 | 10 
Pas [1 | 0 | 3000 | -3000 [18 | 9500 | -72e3 | eer | 441 =| 178 | 12 ~~ 
Paaa tS ~~| 0 =f 2000 | 2000 | 22 | S500 | 74a | 47 | 2100 | 20 | 18 
A A A = 
Pog) 2a OMe OMIT, | 2500, | scam] 01) 02s] 
Paa7 ft 2} Oo | OT CT =O SCL 8 | 9500 | -voms | 01 | 904 | 133 | +16 
Pae ft | (0 =| ~3000 | -s000 | 17 | 9500 | -r2e3 | -eo7 | 4140 | 214 | 15 
Pas [2 | o | -2000 | 2000 | 09 | 9600 | -7eas | 1 | sea | 162 | + 9 
P40, 1 | 7 {| o | o | 31| 9500 | -7203 | soy |  s69 | 374 | 11 
451 
P42 yf 7 | o | 3000 | -3000 | 31 | 9600 | -72a3 | @e7 | 4100 | 276 | 18 
Past 2) ~—(| 1 ~[ ~3000_~+| -3000 [ 14 | 9500 | -7eos | -o1 | s47 | 23 | 9 | 
(45a 5 | 0 | «1000 :~*f 1000 | o7 [ 9500 | -a7vas [| i [ 32200 | so | 5 
pas [5 | 0 | woo | 1000 | 27 | e800 | 74a | a1 [| s214 | 202 | 16 
part 1 [| 71 f{ oO [: 0 |11{ 9500 | -72e3 | eo7 | 683 | 21 | 2 | 
2 | ee oe eo || 500, | scan) eT 900, | 152m | 
Pao tT 5 | o | o | o |29| 9500 | 4a | 41 | 198 | so | 14 | 
Par [2 ~~ ~| o | 1000 | 1000 [| 09 | e500 | -7eas | -o1 | 1948 | 163 | 7 | 
Paes 2 ~—Cd| «CO ~C| «2000 ~*| 2000 | 12 | 9500 | oss | -o1 | sea | i37 | 9 
Pace t 7 ~~Cd| 7 ~«| ~=<3000 +| -3000 {| 09 | 9500 | -72e3 | o7 | 4139 | 353 | 9 | 
Paes [5 CL 1 | «2000 :+| 2000 [| 13 | 9500 | 47a | 41 | 2190 | 369 | 10 | 
2 | = 
ras; 1. | 0 | o | 0 | 08 | 9500 | -7os3 | 67 | 879 | 175 | 17 
ai [71 | © | 1000 | 1000 | 36 | 9500 | -72e3 | -e67 | se7_ | 253 | 15 | 
Pat 2 | 0 | -2000 | 2000 | 38 | 9500 | -7es3 | -91 | 3674 | 242 | 17 i 
Tamayo 2d i «| 0 SC] ~ 0 Cd oe | 00 | ves | or: =~ eSSCOTCCTCOC 
ramets | o | 3000 | -3000 | 16 | 9500 | 478 | 4i [| 5728 | i178 | 6 
crn ee | ON Ome (2500) 728s ce7 | ee | | 285 
Pav [1 {| o | 3000 | -s000 | 22 | 9500 | -72e3 | -e67 | 4738 | 26 | 10 
Tara 5 | 0 | 3000 | -3000 | 09 | e500 | ams | 41 | S73 | o19 | 4 
Tay 2 {| o | o | o | 14 00 | -roa3 | -o1 | e698 | 160 | 11 
ai | 5 | 0 | 9000 | -s000 | 29 | 9500 | #748 | 41 | 5723 | 851 | 11 
[2 LE Ee es ae eee ee ee 
eo ee a Se Se i 
Tas [5 | 7 | 1000 | 1000 | 26 | 9500 | 4748 | 41 | s2i2 | 173 | 16 
Pass [5 | o | 3000 | -s000 [14 | 9500 | ae | a1 | 57 | 550 | 5 
Pee | ng OMI CaM. 7 | S500" (me/cncmlmcoy Imam og) ge eh 
pee | a eee 30 | ae 
489 i i =O | ~ 9500 = | 41 2006 451 
af ne a0} rs tsa 
a 5000 | 000 322500 rasa a sta 
Paez ft 2 | o | 3000 | -a000 | 22 | s00 | -7ea3 | -o1 | 343 | 166 | 14 ‘| 
Pass {2 [| 1 [| 1000 | 1000 [22 | 9600 | -ras | 1 [| ta | 9 | 12 
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case [10] 3000 | -3000 
ef : 
ee ee ce a 
a 
case [ 5-0-2000 | “300008 | 9800} “are [ai] 2108 

499 2 1 2000. | 2000 | 1.7 | 9500 7 SE 
Taso [ 22000 | 2000.7 | as00 | 7689 | at [sera [9 fa 
er eae ee 
300 0005 ss | 588 tf at 
-303-[ 2] 9000 000] 3500} ras | st a 
35 000 F000 8500} res 867 | a 9 
Hs 002000 7 00 8s 867 aie so 
05-000 0008800 | ares at a 
SA I 
Ss a 
i 
cso 2000000 [25 [9500] 7535 a1 | ts a8 
a 
a ti te Hi ee ee ee +B SS 
a 
a == 
Sc ae oe Oe Seo [ares [at [a 

sie | 2 ~~ ~| 1 ~~4{| -2000 | 2000 | 14, ar eo SS 
Sk 
1 
a 
S20 [so 89500] 4788 iss 102 
=e eo gar ee i 
| 522, 1. | 1 ~| 3000 | -3000 | 26 | 9500 | -7283 | -867 | 4140 | i 
ses [aoa a0 ae | 58s] 
524] 2 T1000 -| 1000 | 22 | 9500 | -7503 | -91 | 1347) 274 d= 
sat a1 | 80 a 
z i | “Yo00 | “Yo00 =| 27- [9800 [7a | 367 [sot [a0 a 
A a Se ee ee ee ee ee ee 
car See Raa 24 | 3500 [ras [eer [a8 [af is 
eS 
500 [1 | 0 [3000 | -3000| 4.8] 9500-7283 | 367 | aia 7 
53im| 5 | 1 [mec | ™1Q00 "18 | ©9500 | 4748 1 47 | cone. s+} 
= a OO 
530 7379 
534 [| 1. ~+| 0 | 3000 | -s000 | 31 | 9500 | -7283 | -867 | 4139 ‘| 3 3 
= 
=: Se E 
sar [0 ao00 zc e800 | ras | 387 [ais zis fs 
338 a i 
a : 7 

) a 





Ti ew : 
=f ae ee ik 
MR MM a a 
2 ee eee ee ee 
Ue.) a a a [2000 [2 EO eS ee 
= Z ta 19 9500 -7533 9 1339 ee Ve — 
a 
sar 0 so a8 Se a 
SB 
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Error Control 


Hours from CARP X ; CARP Y ea 
Pred Wind oe ae Predicted CARP | No Control | Controlled | Inputs 


Pee [2 | 0 | Cit [9500-7583 [oi | 903+) Stadt 
ee 
Tesi | 1 | 0 | 3000 |~-3000_| 25 | 9500 | -7283 | e67 | 4143. -[ me | 16 
Tss2 [1.1 | -2000 | 2000 [05 | 9500 | -7283 | -ee7 | 318 | 278 | 14 
sea [1 | 0 | 1000_| 1000 [05 | 9500_| -72es | eer | 473__| 18 | 719 
sea | 5 | 1 | 2000 | 2000 [1.2 | 9500 -a7aa_[—a1_| ato _aao_—sd|~=it 
Dass | 1| 0 | 1000 | 1000 | 04 | 9500 | -72e3 | ee7 | eo1 | 21] 12 
see | 2 [| 0 | 1000 | 1000 [2.4 | 9500 | -7538 [91 | tse7_| 1e4_ 17 
-se7 [5 | 0 | 3000 | 3000 | 1.6 | 9500 | ~4748 | a1 | 5726 | 100d [5 
Tse | 5 | 1 | 3000 | 3000 | 06 | 9500 | 4748 | a1 | s723_| 1304] ~=3 
Tass [2 | 0 | 3000 | 3000 [1.0 | 9500 | -7ea3 | -o1 | 3432 | tei] 11 
seo [2 {1 | 0 | 0 | 09 | 9500 | -75a3 [or | 911 | 252 | 16 
sr; 2 | 0 | 0 | o | 18 | 0500 | 7503] 01 | 009 | 148 | 13 
sez [1 | _1 | s000_| 3000 |-0.6 [9500 | reas | eer | 4142 [tea 3 
ses | 1 |_1 | 1000_| 1000 | 1.2 | 9500 | -72a3 | eer | 586. [| 235_~=~|~C 
[sea 1] 0 | -2000_|~2000 [0.5 | 9500_| 7283 [867 | 3132_| 222_| 10 
ses | 1 | 1_| -2000_| 2000 [2.2 | 9600 | -7283 | s6e7 | 3183__| 27] 3 
ses [1 | 1} 0 | 0 | 24 | 9500_| -7283 | eer | 67 tea_| 11 
Deer | 5} 0 | 0 | 0 | 20 | 9600 |~-a7aa_|_41_|_1007__ | 17__| 13 
see [1 | 0 | -2000 | 2000 [1.7 | 9500_|-7283 | 67 [3136 [ 168 | 12 
see [1 | 0 | 1000 | 1000 | -08 | 9500_| -7283 | -8e7 | S87 203+ 17 __ 
[57 [| 2 | 0 | -2000 | 2000 | 7.9 [9500 | -7ea3 | _-o1 | 3673 | 182) 11 
Csi [| 1_| 1 | -2000-|_2000_[ 1.0 | _9500_|~-7283 | e67 | 3133 270+) 11 
Ps {2 | 1] -2000 | 2000 | 1.6 | 9600 | rss | -o1 | 3673. | 307+) 11 
Tse [5 | 1 | -2000 | 2000 [1.5 | 9600 | -4748_[_41_| 2180 | 265 | 15 
ava [1 | 0 | 1000 | 1000 | 07 | 9600 | reas | 867 | 581 | 29 ~~? 
rss [5 | 0 | 1000 | 1000 [24 | 9600 | 478 | ai | sai7__| 3433 
sve [1 | 0 | 1000 | 1000 [21 | 9500_|_-7263_| 867 [503 | 306] ~—13 
-s77 [2 | 0 | 1000 | 1000_| 1.5 | 9500 | -7sas | 91 | 1954 | 160 | 10 
ree] 1 | 0 | 0 | 0 | 27] 9500 | 7263 | 667 | 806 | 256 [| 13 
reefs | o | o | 0 | 34 | 9500 | 4748 | 41 | 1983 [274 | 20 
Tso [2 | o | o | 0 | o7 | 9500 | ress | or | ase 150 | 11 
Peer [1 | 1 | 3000 | 3000 | 1.8 | 9600 | -7eas | er | 4140_| 216 | 12 
Deez [5 -| 1 | 0 | 0 | 07 | 9500 | 47a _|_ai_| 2001 [525 | 7 
ses [2 | 0 | 1000 | 1000 [26 | 600 [ves | or | 1350 | 173 | 13 
rsea {2 | 1 +{ 0 | 0 | 24| 9500 | -7ex3 | 91 902 | 84] 14 
reese {5 {1 | -2000 | 2000 {| 13 | 9500 | 4748 | 41 | 2103 | 380] 9 
ses [5 | 1 | 3000 | -3000 | 73 | 9600 | -a7aa | a1 | s7aa_ | eaz_| 5 
ears |e? a [eee | OOD [208 A) A om 2500-7 eee Dee 02a © sn el gion 
sf 5 | o | -2000 | 2000 | 28 | 9600 | rae | a1 | 2e7_| 260 | 15 
so [ 1 [| o |] 0 | 0 |16| e600 | -r2e3 | -ee7 [00s | 29 | 10 
RO: a 
Tp S| 1 | 1000 [1000 | 06 [e500 | 474s [41 3zs | ss s=dYC=Cit 
S| sa | a0 22 sano | est 
se [1 | 0 | -2000 { 2000 | 45 [| 0500 | -7263 | 867 | 3130 | 172 | 17 
se] 5 | 0 | woo [i000 [21 | e500 | 47a | a1 | sie7_| 105 | 18 
5] 5 | 1 | 1000 | 1000 | 04 | 9500 | 474s | a1 | s003 | 22 | 1 
pss [1 {0 | 0 | 0 |21 | so | res | ery | ea | 38 | 10 | 
reer [2 | 1 | too [i000 [04 [| e500 | -75a3_| 91) 1352. | 316 | 10 
ss [1 | 0 | 3000 | -s000 | 20 | 0500 | -7283 | 867 | 4139 | 257 | 
ay ed ae 1000 1000 | O02 | 9500 | -7283 | -867 593 18 





i Paes Wid kt Did Dd et DN Fm et i (tn Fie Lh Voha Lad t (hOa Lme )| ink | a) 1d eed PITS IRS ENE ee ean) ay mea sy 


; 600 | a 0 0 foe pee ou as he ayo Flee 1986 217 | 13 | 


On 


103 





SAMPLE SIMULATION RESULTS 


APPENDIX B 


This appendix contains output for numerous simulation tnals to illustrate a sample 


of the results obtained. 


Run35 


Y Error: 


X Error: -2.8414 


-11.3694 Hor Error. 11.7191 


No Control X Error: -466.9893 Y Error. -122.2171 Hor Error: 482.7174 
A Inputs: 4 B Inputs: 3 C Inputs: 6 D Inputs: 2 


Total Control Inputs: 15 


Trajectory Comparison: “Run 3. 


(ee) epninny 


weit 


Distance (feet) 


Bo 
ce 
=: 

= 
So 

oe: 
bis) 
Qa, 
E 

oO 
fee 
So: 
GS 

e 

= 


— Controlled 


rolled 
redicted 


--- Uncont 


-----7-]-—— P 


et ie te ee oe 7 


( 


ia0o}----- 


e9)) 821 


1881p 


2500 


| distance (feet) 


2000 
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ean wale 


a 


~-- Pred Wind N 
Pred Wind E 


=< Uj 
uw 
= 
ss 
SS 
p= = 4 
BB 
Lf 


| 
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Run202 


=2295, 7425 selor Error, (23967913 


X Error: 688.6045 Y Error: 


No Control X Error: 5963.0562 Y Error: -5411.0036 Hor Error: 8052.1425 


A Inputs: 1 B Inputs: 1 C Inputs: 0 D Inputs: 1 


Total Control Inputs: 3 


Ae 
™.: 
ec 
2 
ite 
cS . 
D.s: 
£2. = 
O°: 
Q.: 
&: 
One 
GD eS 


oi 


i Pe fi ee coe ee re oe ee 


™“ 
en 


em ecto m ne te. 


+ 
% 
foav-7= , 
5 


eewewv tf” 


- Distence (ea 


— Controlled 
--- Uncontrolled 


Predicted 


eet cette re arc 


pot ctcecces 


v 


_—ew ew ew ww ow tw] we 


eeeeeeee = 


\ 


we nqnmnnnanen eee ae geen nnee en es 


a 
\ 


< 


wwe eee eo woe owe = 


x 
Y 
‘ 


ve 


weenrtleewea we = 


anaes eae an ook aa aa = 
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Run202 (continued) 


Ee 
ai. 
a 
=. 

am 

ae 


-——— = 


eewrwereew = 
cs 
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Run217 


-1305.1452 Hor Error: 1321.8651 


X Error: -209.5789 Y Error: 


No Control X Error: 7971.5228 Y Error: -7412.1185 Hor Error: 10885.0667 


A Inputs: 6 B Inputs: 5 C Inputs: 5 D Inputs: 6 
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Run217 (continued) 
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Run223 


-37.2997 Hor Error: 121.1473 


No Control X Error: -684.7159 Y Error: 1083.0916 Hor Error: 1281.3755 


A Inputs: 6 B Inputs: 9 C Inputs: 7 D Inputs: 8 


Error -l15.2622-Y Ermor. 
Total Control Inputs: 30 
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Run223 (continued) 


oe 
i: 
a 
ia a 
nae 
a.: 
Pars 
c 
& 
hae 


seu wee 2 @ 


e2sy7 ec @ 
. 


--- Pred Wind N 


Pred Wind E 
—— Actual Wind N 





112 


APPENDIX C. SIMULINK® REALIZATION 


The following presents the Simulink® realization for this effort. It is presented in 
the hierarchy of the implementation of the model. 


MUST LOAD pmdata.mat BEFORE RUNNING 
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APPENDIX D. AGAS INSTRUMENTATION DESCRIPTION 


Overview 


The AGAS instrumentation package was required to 
measure and record the state of the AGAS payload to 
include: position, velocity, acceleration, and three- 
axis attitude and attitude rates. In addition, the state 
(extended or retracted) of each of the four actuators 
was recorded. The position and velocity of the AGAS 
payload was provided by a differential carrier-phase 
GPS system known as the Improved Vehicle 
Tracking System (IVTS). Acceleration data were 
derived from a tmad of accelerometers, and the 
attitude of the package was measured by an Attitude 
Heading Reference System (AHRS). Pressure 
transducers attached to the PMAs provided the state 
of the actuators. Since time correlation of the data 
from multiple sources was critical, a timing source 
was included in the package. All measured data were 
recorded on removable flash memory devices. 


The AGAS instrumentation package contained two 
PC/104 computer systems. The sensor computer 
system managed the data acquisition and recording 
for the AHRS, acceleration and actuator states 
collected at a 15-Hertz (Hz) rate. The Global 
Positioning System (GPS) computer system 
controlled the GPS receiver and acquired and 
recorded the receiver measurements. The GPS 
trajectory solution was recorded at a 5-Hz sate. 
Details of the instrumentation system design are 
presented in the appendix. The system design is 
illustrated in Figure 1. 


Sensor Instrumentation GPS Instrumentation 





Figure 1. Instrumentation Overview 


System Architecture 


The data acquisition system was based on an industry 
standard embedded computer product family called 
PC/104. This product family consists of the central 
processing units, power supplies, and video as well as 


a large assortment of analog and digital data 
acquisition and contro] modules useful in configuring 
data acquisition systems. The CPU modules are 
based on Intel® X86 processors and use the Industry 
Standard Architecture (ISA) bus associated with the 
IBM® PC/AT series of personal computers. Typical 
enhancements to the standard computer include solid 
State mass storage, Basic Input/Output Systems 
(BIOS) extensions, and reduced power requirements. 
These modules are highly integrated and have a small 
footprint, measuring 3.6 X 3.8 inches. The stacking 
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Figure 2. PC-104 Modules 


feature of the PC/104 family makes it easy to 
configure a system to meet unique requirements. A 
typical PC/104 system is shown in Figure 2. 


AGAS Instrumentation Computer System 


The sensor computer system contained a central 
processing unit (CPU) module, an analog-to-digital 
converter (ADC) module, a power supply module, a 
timing module and a type II PCMCIA module. The 
CPU module, is a highly integrated module 
containing an Intel 486 DX4 processor operating at 
100 Mhz, an IDC hard disk controller, a floppy disk 
controller, two serial ports, a parallel port and 4 
megabytes of random access memory. The module 
also contained a 1 megabyte flash memory which, 
through BIOS extensions, looked to the processor 
like a bootable IDC hard disk. This flash memory 
was programmed with the operating system and the 
application software that controlled the modules on 
the stack. 


The ADC module is an 8 channel 12-bit converter 
configured for 5-volt bipolar input An internal 
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adjustable gain amplifier provides programmable | 


gains of 1, 10, 100, or 1000. The analog to digital 
conversion rate is a function of the mode used to 
initiate the conversion, and varies from 2,000 
conversions per second in the software conversion 
mode to a maximum of 100,000 conversions per 
second if the direct memory access (DMA) 
conversion mode is used. An interrupt conversion 
mode is also available, which will provide up to 
20,000 conversions per second. In this mode, the 
interrupt signal may be provided by an external 
signal or it may be provided from the output of a 
programmable 32-bit counter contained on the 
module. This module also contains two 12-bit 
digital-to-analog converters as well as four digital 
inputs and four digital outputs. 


The timing module is a digital clock that provides 
accurate time to application software through the 
PC/104 bus. The clock may also be initialized from 
the application software through the bus. Time-of- 
year is available in a binary coded decimal (BCD) 
format with a resolution of hundreds of nanoseconds. 
Clock accuracy is maintained by phase-locking the 
internal clock to an external reference. The reference 
mav be any of the standard timing signals usually 
found on United States test ranges, such as the Inter- 
Range Instrumentation Group (IRIG) time code 
signals or NASA time code signals. In addition the 
internal clock may be synchronized to an external 
one pulse per second (1-PPS). Synchronization 
accuracy to the IRIG or NASA timing signals is 
specified to be between 20 microseconds and 5 
microseconds depending upon which code is being 
used for the synchronization input. A 
synchronization accuracy of one microsecond is 
specified for the 1-PPS input. 


The sensor computer system power supply module is 
a DC to DC converter that provides up to 50 watts of 
power in the form of +5 volts and +- 12 volts. The 
power supply module requires a DC input voltage 
between 6 volts and 40 volts. Each of the power 
signals is connected to the PC/104 bus to facilitate 
powering the modules in the stack. 


The PCMCIA module located in the sensor computer 
module stack is a standard type II PCMCIA interface. 
With the appropriate driver software, flash memory 
modules inserted into the interface appear to the 
operating system as an IDC hard drive. 


Three modules compnsed the GPS system computer: 
a CPU module, a power supply module, and a 
PCMCIA module. The hardware on these three 
modules is identical to the hardware in the sensor 


system computer modules. The CPU was 
programmed with the software required to interface 
with the GPS receiver sensor. A photograph of the 
packaged instrumentation system is included in 


Figure 3. 
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Figure 3. AGAS Instrumentation Package 
AGAS Sensors 


As mentioned, the AGAS instrumentation package 
contained four types of sensors: a GPS receiver, an 
AHRS unit, an accelerometer triad, and muscle state 
sensors. The GPS sensor consists of an L] C/A code 
receiver and antenna. The receiver contains 12 
independent channels, which track the GPS satellite 
signal in parallel. The receiver is capable of making 
measurements at a 20 Hz rate and providing position 
and velocity output up to a 10 Hz rate (a 5 Hz rate 
was used for the AGAS implementation). Two RS- 
232C VO ports are available on the receiver to 
provide for control inputs and data outputs. Either of 
the I/O ports may be used to input commands and 
output data, or both may be used at the same time. An 
auxiliary I/O port provides a 1PPS. 


The AHRS used in the AGAS instrumentation system 
is a solid state device that measures three-axis 
attitude rate and integrates it to form attitude and 
heading. The attitude and heading is compared with 
two reference vertical pendulums and a tri-axial flux- 
gate magnetometer. The resulting error is filtered 
and used to adjust the output of the svstem causing it 
to converge to the attitude of the vertical pendulums 
and to magnetic heading. The attitude and attitude 
rates are available on an RS-232C interface at a rate 
of just less than 15 Hz. The specified accuracy of the 
AHRS data is shown in Table 1. 
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Dynamic 
Static | (percent) 
Rate Accuracy (deg/sec) 


Anitude Accuracy (deg) 
Heading Accuracy (deg) 


Table 1. AHRS Accuracies 

The accelerometers in the AGAS instrumentation 
system are miniature devices with a full-scale output 
of + 10g. The acceleration experienced by these 
devices is converted to an electrical signal via strain 
gauges wired as a Wheatstone bridge and attached to 
the acceleration-sensing element. The frequency 
response of these devices is nominally 0 to 300 Hz, 
with a minimum response of 140 Hz. The sensitivity 
of each element is approximately 8 milli-volts/g. 
Three accelerometers were mounted in a steel 
mounting block with their input axis orthogonal. One 
axis was defined to be the X-axis with the + input of 
the accelerometer in the direction of +X. This axis 
was designated as the axis pointing forward. The 
other two devices were mounted such that the +Y 
axis was pointing to the nght and the +Z axis was 
pointing down. The mounting block was attached to 
the mounting structure with the +X axis of the 
accelerometer triad parallel to the AHRS +X axis. 





The last group of sensors in the AGAS sensor system 
is the riser muscle state sensors. These were simply 
four pressure transducers, one for each muscle. The 
four muscles in the control system have two states, 
extended and contracted. When the muscle is 
contracted, the pressure in the lines supplying air to 
the muscles is about 50 psi. When the muscle is 
extended, the pressure drops to nearly zero. 
Therefore, the pressure in the lines represents the 
state of the muscles (zero indicating the control is 
activated). 


AGAS Instrumentation Software 


The software residing in the processors shown in the 
diagram controlled the hardware in _ the 
instrumentation package. Both processors ran the 
MSDOS® operating system, and the application 
software was written in the C language. 


The application software in the GPS processor was 
the same as that used in the Yuma Proving Ground’s 
Improved Vehicle Tracking System, and in 
conjunction with the receiver used, will produce sub- 
meter positioning when the data is post-processed to 
remove errors. This software initialized the GPS 
receiver at power on and began to acquire the GPS 
satellite signal. While the AGAS package was in the 
aircraft, the GPS receiver tracked the signal via a 


retransmission system that received the signal 
external to the aircraft and retransmitted it inside the 
aircraft. Once the GPS receiver was tracking 
satellites, the processor collected the navigation data 
and recorded it on the recorder. 


The sensor instrumentation processor also began 
initializing at power up. Early in this process, it 
started requesting a time message from the GPS 
receiver. Once this message was received with status 
indicating the GPS receiver had solved for GPS time, 
the sensor instrumentation processor initialized the 
timing module to the current UTC time. The 
accuracy of the time in the timing module was 
subsequently maintained by phase locking its clock to 
a 1 PPS signal from the GPS receiver. Once accurate 
time had been established in the timing module, the 
sensor processor enabled the AHRS output and began 
acquiring data. 


The attitude and attitude rate data from the AHRS 
was transferred to the processor through an RS-232C 
interface. Each transfer contained a string of ASCII 
characters that included a unique header, the current 
AHBRS attitude and attitude rate. The transfer was 
terminated with the ASCII character representing a 
carriage return. A period of inactivitv on the 
interface separated the transfers. At the point the 
software detected the unique character designating 
the start of a transfer, the current time was sampled 
and stored. This time became the time of the attitude 
sample. 


As soon as the attitude time sample was stored. the 
software acquired the accelerometer and pressure 
data from the ADC. This activity started with 
another sample of current time, which became the 
sample times of the data from the A/D converter. 
After the accelerometer and pressure data were 
sampled, the software converted the data to ASCII 
and recorded it in a file designated for A/D converter 
data. 


At this point the software continued with the 
acquisition of the attitude record. Since the serial 
port driver: operates under interrupt control, the 
attitude data that arrived at the port during the A/D 
sample interval was placed in a buffer. These data 
were retrieved, formatted, and recorded in the file 
designated for attitude data. This finished one loop 
of the data acquisition process and the software 
began to monitor the serial port for the beginning of a 
new attitude sample. The message rate from the 
AHRS was approximately 15 Hz, so the sample rate 
for all data was also at approximately 15 Hz. 
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APPENDIX E. MATLAB” SCRIPT FILES 


The following presents the MATLAB® script files utilized in this effort. Note 


that only the critical software code is presented. 


Initialization File 


$This program initializes the workspace for running the 
simulation. 


load pmdata 


EUNnS-1nput {(*Desired Number of Runs: ‘'); 


disp('To accept current parameters (one by one) just hit 
Eccurn.);adisp(* “);displ(* “) 


disp('Current Windfile: ');disp(windfile) ; 
newwindfile=input ('Select a Wind File ([day drop#] eg. [310 1]}): 


ye 


newwindfile 


if isempty (newwindfile) ~=1 
windfile=newwindfile; 


end 
for i=l:runs 


if windfile (1,1) ==310; 
RAND=ceil (rand(1)*3.999); 
filechoice=[110607 110609 110612 110613]; 
end 
if windfile(1,1)==311; 
RAND=ceil (rand(1)*2.999); 
filechoice=[{110707 110711 110714); 
end 
if windfile(1,1)==0617; 
RAND=ceil (rand(1)*3.999); 
filechoice=[061710 061712 061714 061716]; 


end 


12] 


rawinfile=filechoice (RAND) ; 
saon—t loor (rand(1)*1.999) ; 
$sastart=ceil (rand(1)*6.999) ; 
RAND=ceil (rand(1)*3.999); 
GtEsece=[0 0; TUCO" 1000 ; -20009 2000; 3000 F-s0c0. 
xoffset=offset (RAND,1) ;yoffset=offset (RAND,2); 
psidot=randn (1)+1.85; 
BracDir=t loom (land ai 22): 
tee eash ie — 
CompassSign=1; 
else 
CompassSign=-1; 
end 


CompassBias= (rand(1)*1.999) *CompassSign; 


Le oeeeteba ah at aeee 1) as any 
releasealt=20000; 

else 
releasealt=9500; 


end 


te eextist ( amhirs le') 0 
initfile=[windfile rawinfile saon releasealt Vi Hi dzalt 
xoffset yoffset psidot CompassBias]; 
Save init initfile 
else 
initfile(i,:)=[windfile rawinfile saon releasealt Vi Hi 
dzalt xoffset yoffset psidot CompassBias]; 
save init initfile -append 


end 
clear new* acchange 


end 


disp('init done') 


2 


Equations of Motion — dof3.m 


function accel=dof3(z); 
V=2(1:3); 
Vt=norm(V); 
%checks to see if density is contained in the input vector. 
%If not, sets to standard sea level density 
if length(z)==5; 
rho=2(4); 
time=z(5); 
else 
rho=0.002377; 
time=z(5); 
end 
%Atmospherics 
qbar=.5*rho*(Vt*2); 
g=32.17; 
%System aerodynamics/characteristics 
opentime=3; %seconds 
Gd=0:75: %obtained from flight test 
Do=28; %reference diameter in feet 
Dp=.67*Do; %profile diameter 
So=pi/4*Do%2; %reference area (flat circular parachute) 
W=346; 
m=W/g; 
“Calculate drag area. Assumes a linear opening of the parachute 
if time<=opentime 
CdSt=[0 0O;opentime Cd*So]; 
CdSo=interp1(CdSt(:,1),CdSt(:,2),time); 
else 
CdSo=Cd*So; 
end 


“Calculate mass terms 
alpha11=0.25*rho*4/3*pi*(Dp/2)%3; 
alpha22=alpha11; 
alpha33=2*alpha1l1, 


Mi=[(m+alpha11) 0 0; 
0 (m+alpha22) 0; 
0 (m+alpha33) J; 


saccel=[-qbar*CdSo/m*V (1) /Vt;-qbar*CdSo/m*V (2) /Vt; - 
gqbar*CdSo/m*Vv (3) /Vt+g] ; 


accel=inv (M1) * ([-qbar*CdSo*V (1) /Vt; -qbar*CdSo*V (2) /Vt >... 
-qbar*CdSo*V(3)/Vt]+[0;0;W]); 





Monte-Carlo Simulation carlo.m 


%This file executes the desired simulations 
clear 
load pmdata;load init;[row col]=size(initfile); 


Startrun=input('Key-in Starting Run Number '); 
for runno=1:row 
windfile=initfile(runno,1:2); 
eval(['Wind=Wind' num2str(windfile(1)) 'd' num2str(windfile(2)) ';']) 


rawinfile=initfile(runno,3); 


eval([‘rawin=rawin’ num2str(rawinfile) ';']) 


saon=initfile(runno,4); 
if Saon== 
saseed=[floor(rand(1,3)*1000) zeros(1,6)]; 
else 
saseed=floor(rand(1,9)*1500); 
end 
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releasealt=initfile(runno,5); 

Vi=initfile(runno,6); 

Hi=initfile(runno, 7); 

dzalt=initfile(runno,8); 

xoffset=initfile(runno,9); 

yoffset=initfile(runno, 10); 

psidot=initfile(runno,11); 

CompassBias=initfile(runno,12); 

sim(‘carp’') 

disp([num2str(runno) ' CARP Sim(s) Done'}); 
CARP=-carptraj(length(carptraj),1:2); 
fliptraj | 

sim('c9pointmass') 

%noisel =wavread('c:\windows\media\office97\gunshot.wav'); 

%sound(noise1) 

disp({[num2str(runno) ' Control Sim(s) Done'}) 

monte; 

clear carptraj; 

horimpacterr 

results(runno,:)=[runno+startrun-1 windfile rawinfile saon ... 

xoffset yoffset psidot releasealt Vi Hi dzalt CARP ... 


hornocontimpacterr horimpacterr totcontrols]; 


end 


eval({'save runfile’ num2str(startrun) '_' num2str(startrun+runno-1)... 


"asc results -ascii'}) 
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